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Introduction 


This final report suimnsrises the results of the analysis of data acquired 
by two flights of an array of six Bragg crystal spectrometers. The objective 
of the experiments were to obtain high spatial and spectral resolution observa- 
tions of various coronal features at soft X-ray wavelengths (9-23^). 

The first rocket, NASA-Aerobee 26.002, was launched on 11 June, 1973, 
at 18:30 U.T. It observed two active regions and an area of quiet corona. 

Both active regions, McMath 12379 and McMath 12387, were quiescent at the 
time of observation. The region of quiet corona was near the center of the 
visible disk. The second flight, NASA-Aerobee 26.003, was launched on 19 
December, 1973, at 22:07 U.T. This observed a single active region, McMath 
12664. This region was fairly active prior to the flight. 

Each spectrometer system consisted of a two-dimensional collimator with 
a field of view of approximately 35 arc seconds full width half maximum 
(FWHli). Two of the instruments had larger fields of view of 150 arc seconds. 
The diffraction crystals were KAP, and were mounted on platforms that were 
rocked through the desired angular ranges by a motor driven cam. The position 
of the platform was monitored by a shaft encoder that divided the scan into 
512 angular bins. The detectors were gas filled proportional counters. The 
X-ray events were summed in a register, and along with the various other 
instrument functions were recorded in the ground receiving station via the 
PCM telemetry. This was repeated at a rate of 256 Hz. An H-alpha camera re- 
layed pictures to the experimenters who were able to control the aspect of the 
rocket. The images were also recorded on a 16mm film camera for post flight 
analysis. For further details see Final Report, Contract NAS2-6723 dated 
15 February 1976. 
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In the report we describe the various aspects of the analysis of the X-ray 
data obtained by the two flights. These observations were coordinated with 
observations from the experiments on the Apollo Telescope Mount (ATM) . Part 
of the work done during the course of this investigation was to relate the 
various data sets to one another. We participated in the Skylab Solar Workshop 
Series on Active Regions in Boulder, Colorado, which was sponsored by NASA 
and NSF. The Apendices attached to this report contain the published results, 
abstracts of papers, computer code descriptions and preprints of papers, 
all produced as a result of this research project. 
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Scientific Results 


(a) Background 

The basic goal of this research project was to obtain a holistic model of 
a solar active region using data from a wide range of instruments. The 
opportunity arose to tackle this problem with the coincidence of the flights 
of the two LPARL crystal spectrometer payloads with active observation periods 
on Skylab’s ATM. 

The project naturally divided itself into three distinct phases. The 
first consisted of the analysis and interpretation of the soft X-ray emission 
line data from the rocket borne experiments. The reduction of these data 
was essentially complete prior to the start of the contract. However, the 
interpretation and analysis of these data brought to light some interesting 
results. The second phase was to correlate the rocket observations with the 
cotemporal observations from other in “ *’rum^nts , primarily those on the ATM. 

This was the most difficult part of the project to complete as it involved 
codifying the observations from other investigating groups, defining and cal- 
culating a common aspect solution as well as determining a relative calibra- 
tion. This stage still remains the single largest possible source of un- 
certainty in the analysis. The last phase of the project involved modelling 
the regions using the assemblcJ data and comparing the results. These 
stages are described in the following three subsections. The concluding 
subsection discusses briefly some of the interesting ,iceas that have resulted 
from this project, and especially during the Skylab Workshop on active regions, 
and the primary lessons learned from the analysis. 
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( b ) A nalysis of the X-Ray Emission Line Data from the Rocket Borne Experiments, 
A vast quantity of data was accumulated during the two flights of this 
experimental package. Spectra were recorded from 31 stable pointing positions 
by all six spectrometers. At each pointing position every spectrometer pro- 
duced two spectra of its dedicated wavelength band, one scanned in the wave- 
length increasing direction and the other scanned in the reverse direction. 

A total of 58 X-ray emission lines were definitely resolved with a possible 
20 or more marginal features. Figure 1 shows typic ;1 spectra from the 
wavelength bands monitored during the flights. 

The reduction and analysis of these data are summarized along with the 
results and will be submitted for publication in Solar Physics by Strong 
and Acton (Appendix 1). In Appendix 2 a full set of the results of the 
analysis are listed. As a result of the analysis we are able to summarize ’ 
the primary results and conclusions under the following subheadings: 

(i) Spectral Line Fitting Techniques. 

We have developed a technique that fits an entire spectrum with a 
series of line profiles superimposed on a constant background using an 
iterative least squares technique. A Lorentzian line profile was assumed 
to represent our data. In the fitting procedure the position, width and 
peak count rate of each line as well as the background were left as free 
parameters where possible. Figure 2 shows the fit of our Fe XVII cpectral 
interval with such a function. 
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This technique made it possible to derive the integrated line 
intensity for each line from the fitted function. This has several advantages 
over summing the data and subtracting the average background which is an often- 
used method. Firstly, it eliminates confusion between any contribution to the 
background from the continuum and the wings of the lines in the spectrum^ 

Secondly, by using the widths and positions of the more prominent lines it 

is possible to derive relatively accurate Integrated line count rates for faint 
or blended lines. The uncertainties on the characteristic parameters of a 
given line derived from the least squares fitting procedure makes it possible 
to estimate quantitatively the uncertainties in the line count rates. 

(ii) Preflight Calibrations. 

We had available a comprehensive set of preflight calibrations 
which Included the transmissions of the collimators, window support grids and 
detector windows. The integrated reflectivity of the crystals were kindly 

measured for us by A. Bureck. By combining these efficiencies, it was 

possible to derive the integrated line intensity for each line with the 
associated uncertainty. 

(iii) Line Ratioing 

Understanding the sources of the uncertainties on the line inten- 
sity made it possible to minimize their effect on any given result. The errors 
will be of two basic types. There are random errors introduced by statistical 
effects on measured quantities as well as systematic errors resulting from 
instrumental effects (geometry, wavelength, etc.). Also, there are problems asso- 
ciated with atomic data which are not always reliable. The best way of minimizing the 
selective effects is to take ratios of carefully chosen lines where some of 
these factors may be eliminated. The ratio of two emission lines for which 
emission can be represented by an effective temperature and emission measure is:- 
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adopting the usual nomenclature. Chosing lines from the same spectral scan can 
eliminate relative instrumental calibration error and reduce wavelength 
dependent effects. Other factors can be eliminated by choosing the correct 
combination of lines, for example:- 

Two lines of the same ion yield a ratio dependent on their 
relative oscillator strength 




5040 (W^-W^) 


'eff 


( 2 ) 


Ratioing tvo lines of adjacent ionization stages of the same element will 
give information on the equilibrium state of the emitting plasma and its 
tempe rature . 
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If we choose two lines with similar temperature functions then 


5040 (W-WJ 

N -( j- -— ) 

C(T) • TT" ’ 10 eff ^ constant, 
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then we obtain a ratio that is sensitive to the relative elemental abundances 



Three typical examples are calculated in Tables I - III. 
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These methods were applied to our data with some useful results. Especially 
interesting were the results of the abundance ratio of OrN^. We had sufficient 
results to do a statistical study of the problem. The data basically gave a 
result of 6.0(0. 5) however, there were some interesitng departures from the 
value. Some regions from the flight of Aerobee 26.002 indicated a value of 
12(2) assuming active region temperatures ('^ 3.10^ °K) . These results can 
be explained in terms of a higher temperature component in the emitting plasma 
with some justification from the Ha observations from Big Bear Observatory. 
However, the strongest region observed in McNath 12387 gave a value of 4. 0(1.0) 
which can be explained in terms of a real difference in the abundance 
ratio if the observations are correct. Figure 3 shows the results of this 
particular study. See .appendix 1 for a more detailed discussion of this analysis. 

(iv) Resonance Scattering 

Probably one of the most important results of this analysis has 
been the discovery of significant variations in the ratios of the sum of the 
forbidden and intersystem line intensities to that of the resonance line 
(the G ratio) for the 0 VII ion (see Figure 4) . We have been able to interpret 
these data in terms of resonance scattering. This has triggered a study of the 
problem comparing a simple model of the scattering of X-rays in the corona 
with observations of the equivalent case. The resulting paper has been sub- 
mitted to the Astrophysical Journal, a copy is contained in Appendix 3. We 
have also put a proposal to NASA for an expanded study of the problem and its 
implications. 

(v) Modelling of the Active Regions. 

The relative brightnesses of the most prominent lines from each 
stable pointing position in and around each active region made it possible to 
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TABLE I 


Ttenperature Dependence of the 0 VIII (18.95A) to Ne IX (13.44A) 

Intensity Ratio (Jordan) 


^6 

Log T 

1.352 

■^6 

F, > 

1.00 

6.00 

1.352 

-0.085 

1.26 

6.1 

1.073 

-0.085 

1.58 

6.2 

0.856 

-0.085 

2.00 

6.3 

0.676 

-0.085 

2.51 

6.4 

0.539 

-0.085 

3.16 

6.5 

0.428 

-0.085 

3.98 

6.6 

0.340 

-0.085 

5.01 

6.7 

0.270 

-0.085 

6.31 

6.8 

0.214 

-0.085 

7.94 

6.9 

0.170 

-0.085 

10.00 

7.0 

0.135 

-0.085 


^^^Je IX^ 

VIII ^ 

Ratio 

-0.03 

-2.09 

0.161 

-0.01 

-1.28 

0.522 

0 

-0.72 

1.124 

-0.01 

-0.42 

1.517 

-0.04 

-0.36 

1.361 

-0.15 

-0.53 

0.918 

-0.38 

-0.82 

0.653 

-0.80 

-1.14 

0.700 

-1.33 

-1.46 

0.999 

-1.90 

-1.75 

1.719 

-2.44 

-2.00 

3.093 
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TABLE II 


Ttenperature Dependence of the Ne X ( 12 . 138 ) to XI i9.17) 

Intensity Ratio 


T 

Log(T) 

1.667 

T 

1 1 
Log(=i • 3 ^) 

2 ^2 


xr^ 

Ratio 

1.00 

6.0 

1.67 

-0.130 

-5.16 

-1.09 

0.003 

1.26 

6.1 

1.32 

-0.130 

-3.75 

-0.30 

0.005 

1.58 

6.2 

1.05 

-0.130 

-2.63 

-0.07 

0.023 

2.00 

6.3 

0.84 

-0.130 

-1.71 

-0.02 

0.105 

2.51 

6.4 

0.66 

-0.130 

-1.02 

-0.01 

0.331 

3.16 

6.5 

0.53 

-0.130 

-0.58 

-0.01 

0.676 

3.98 

6.6 

0.42 

-0.130 

-0.36 

• 

0 

1 

0.933 

5.01 

6.7 

0.33 

-0.130 

-0.40 

-0.11 

0.813 

6.31 

6.8 

0.26 

-0.130 

-0.62 

-0.24 

631 

7.94 

6.9 

0.21 

-0.130 

-0.92 

“0.60 

0.575 

10.00 

7.0 

0.17 

-0.130 

-1.20 

-1.06 

0.794 

12.59 

7.1 

0.13 

-0.130 

-1.49 

-1.58 

1.230 

15.84 

7.2 

0.11 

-0.130 

-1.74 

-2.19 

2.692 
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TABLE III 


Temperature Dependenae of the Mg XI (9.17) to Fe XVII (15.01) 

Intensity Ratio 


Log(T) 


-2.651 






Xl’ ^<*FeXVIl' 


1.0 

6.0 

-2.651 

1.26 

6.1 

-2.104 

1.58 

6.2 

-1.678 

2.00 

6.3 

-1.326 

2.51 

6.4 

-1.056 

3.16 

6.5 

-0.839 

3.98 

6.6 

-0.666 

5.01 

6.7 

-0.529 

6.31 

6.8 

-0.420 

7.14 

6.9 

-0.334 

10.00 

7.0 

-0.265 


-0.681 

-0.681 

-0.681 

-0.681 


-0.681 

-0.681 

-0.681 

-0.681 

-0.681 

-0.681 


-1.09 
-0.30 
- 0.02 
- 0.01 
■ 0.01 
-0.04 
- 0.11 
-0.29 
-0.60 
-1.06 
-1.58 


-4.77 

- 2.00 

-0.44 

- 0.10 

-0.04 

-0.17 

-0.37 

“0.76 

-1.41 

-2.36 


48.32 

0.418 

0.026 

0.023 

0.034 

0.052 

0.074 

0.114 

0.216 

0.682 


-0.681 


create simple maps of the emission using the higher spatial resolution data 
35**), see Figure 5. Generally we found a close correspondence between the 
resulting maps for various lines with a tendency for the lower temperature lines 
(0 VII and 0 VIII) to have a t:jore extensive area of emission. 

Using a modelling technique first described by Withbroe (1975) 
we derived a series of models of the differential emission measure for each 
stable pointing position. Generally the X-ray emission line intensities could 
be reproduced by an isothermal peak (2.0 - 4.0.10^ ^K) with a width of 
0.3 - 0.4.10^ ^K. This result could be expected as the thermal conductivity 
of the coronal plasma is high at these temperatures, and the temperature 
resolution of the X-ray line is low. 

Th.e X-ray lines available to us from these spectra are only really 
sensitive to temperature in the range 1.5.10^ (0 VII) to 8.0 10^ ^K. 

(Ne X). Hence, we have little or no information on the cooler plasma 
associated with the transition region or particularly hot material from flaring 
regions. Hence, there was a need to call on other sources of data to fill 
these gaps in our models. 

We did carry out an investigation of the stability and applicability 
of the Withbroe technique in the light of some criticism of Craig and Brown 
(1976, and subsequent papers). We found the results to be generally reliable 
and not to reflect the condition numbers of about 50 indicated by Craig and 
Brown. 
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(c) Modelling of McMath 12379 and 12387 


Sinultaneous observations of McMath 12379 and McMath 12387 were available 
from the instruments on the Apollo Telescope Mount (ATM). Hence there were 
available to us a wide range of X-ray and EU\' data on these regions. These 
instruments offered an expanded temperature range for the models into the 
lower corona ( ;^3.10^ °K) as well as an improvement in the temperature resolu- 
tion of the model due to the more narrow temperature functions of the EUV 
lines . 

There were four ATM experiments used in this investigation. S0-5A (AS&E) 
took high spatial resolution (y 2") filtergrams in the soft X-ray (2-5oX) part 
of the solar spectrum. SO-56 (MSEC) was a similar instrument. These were 
used to reconsile the aspect solution of the X-ray data from the rocket borne 
spectrometers as well as be a check on the final models derived from the X-ray 
and EUV emission line data. S055 (HCO) produced raster images (5 x 5 ) at 
a spatial resolution of 3" in several prominent EU\^ lines. S0-82A (NHL) 
produced overlapping images of the solar disk at a 1" resolution in the 180 - 
65oS wavelength band. Data from these experiments were used in conjunction 
with the X-ray line intensities from the rocket flights to produce the dif- 
ferential emission measure models. 

It is evident from the descriptions of these instruments (for greater 
details on them see "Applied Optics" special edition on Skylab) that there 
is a considerable variation in the type of data obtained from the ATM instru- 
ments. Consequently, each had to be reduced, analyzed and correlated with 
our data in different ways. The basic problems can be summarized as follows: 




U) Timing of Observations; not all of the experiments took data at 
the same moment as the rocket data was recorded. Also the observations took 
different times to be accumulated. Consequently, we had to study the preceding 
and subsequent data on these regions to establish whether they were stable 
over the observation periods. However, where possible, data not taken during 
the time frame of the rocket flights were not used. This is especially important 
for McHath 12664 which was particularly active. 

(ii) Spatial Resolution; each data set had its own definable spatial 
resolution. This varied from about 1 arc second (NRL-S082A) to the rocket 
data ( 35 arc seconds FWHM) . Hence, the data sets had to be made compatible 
with one another. As we can do very little to improve the resolution of the 
X-ray data taken through the collimators on the Bragg crystal sepctrometers , 
then it was necessary to degrade the spatial resolution of the ATM data to 
match it. This was done by multiplying the collimator efficiency with the 
intensity of radiation from each resolution element of the image as a function 
of position from the calculated point location. 

(iii) Aspect Solutions; we had to relate the calculated locations of 
each stable rocket pointing position to a position on the image raster of each 
experiment. The absolute pointing parameters for ATM were not well established, 
nor were their relative pointing at least to tht accuracy required in this 
investigation. 

(iv) Calibration; we had to establish the relative calibrations of the 
various ATM instruments and the X-ray crystal spectrometers. 

The ATM instruments were not designed with this type of analysis in mind. 
There have been attempts to combine X-ray emission line data with X-ray filter- 
grams (Davis et al., 1975; Pye et al., 1977). While the temperatures of the 


regions studied have been in good agreement the total emission measure has 
shown some discrepancies. These investigations Involved the use of data 
from a wide field of view spectrometer (3-4 arc min. FWHM) . The finer 
resolution available to us made the resolution of the aspect problem vitally 
Important to this investigation. The analysis is summarized under the following 
subheadings: 

(i) The Available Data. 

Table IV lists all the sources of data available to this investi- 
gation for both flights. 

(ii) Solving the Aspect Problem. 

As mentioned above, this problem was the most important (and time 
consuming) in being able to compare these data sets. The rocket aspect 
solutions were derived by Acton using the outputs of the Fine Sun Sensor (FSS) 
and correlating them with the photographic data recorded by the on-board 
H-alpha camera. We believe these to have an uncertainty of ^ 20 arc seconds 
in absolute terms. However, the relative pointing positions (i.e. the pattern 
of the stable pointing positions) over each region is more accurate. Figures 
6(a) and (b) show the final rocket aspect solutions with respect to the solar 
disk and active regions. 

We used the SO-54 images to establish pointing locations of the rocket data in 
conjunction with the relative intensities of the Fe XVII emission line at 15.01A to 
improve or, at least, confirm our pointing with respect to the ATM instruments. We 
noted that the Fe XV and Fe XVI images from the SO-82A data were similar to those 
images of the "softest" filter available to SO-54 (Filter 2 - patrol normal) . 

Using the IDAPS system at MSFC a pointing position pattern of collimator 
functions were moved over the regions observed until a best fit was found 
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TABLE IV 


Sumrary of Coordinated Cfcservations Associated with LPARL Rocket. 


Observations ■ 

26.002 

26.003 

■ ATM Rocket launch Time 

18:30 UT, 11th June 1973 

22:05 OT, 19th December 1973 

1. ATM Instruments 



AS&E spectrographic 

Good coverage of both 

Thick filter only. 

Telescope S054 (4-60A) 

AR's 

EUV spectrcneter-spectrO“* 

12387 only, grating 

Good coverage 

helicneter SO-55 
(280-1350A) 

in wrong position 


X-Ray Telescope S056 

Good coverage of both 

Good cov^erage 

(3-53A) 

AR's 


XUV spectroheliograph 

Good coverage of both 

Prior coverage may be 

S082A (171-630A) 

AR's 

useful 

UV Spectrograph S082B 
(970-1970A) 

Good coverage of 12387 

7 

White light Coronagraph 
SC^52 


7 

Ha Telescope (6563A) 

Good coverage 

Good coverage 

2 . Ground Based 



Big Bear (Ha Full Disk) 

Sac Peak (Full disk 
spectrcheliograms) 

18:16 - 02:15 UT 

16:27 - 23:28 

D3, H 0 - 4 . 3 IA & Ca k 

18.09 

21.07 

Ha, 3 + Ca k + 0.38A 

18.12 

20.55 

Ha + 0.5A 

18.18 

21.02 

Big Bear (McMath 12387) 



Ha-0.5A, 4. 2x5.6 arc min 

18.12 - 02.05 

X 

Ha-0.5A, 4. 2x5.6 arc min 

18.12 - 02.15 

16.27 - 23.28 

CaK-0.5A, 4. 2x5.6 arc min 

18.12 - 02.05 

X 

CaK-0.5, 1/6 Disk 

X 

17.50 - 22.40 

Lockheed Filter grams - 



large scale 



He I D3(5876A)-8 x 6 
arc min 

18.26 - 18.41 

X 

Ha Doppler (Ha+0.6A) 

10x7.5 1 

17.33 - 22.01 

X 

HA Multislit spectro- 

17.38 - 22.01 

X 

graph 10x7.5 

< 
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Cbservatlons 

26.002 

26.003 

Nbgnetograms 



Crimea 

McMath 12379 

X 

Coronographs & 
Proranences 
Maung Lou (6200A) 

17.00 - 03.00 

17.00 - 03.00 

Sac Peak 

18.43 - 20.22 

X 

Radio Observations 
Clark lake 20-120 MHz 

15.18 - 02.18 

X 

15-2. 5m 

John Hopkins Applied 
Physics Lab 565 1000 MHz 

12.19 - 19.56 

X 

Harvard 10-4000 MHz 

X 

1 

sunrise - sunset 



between the normsllzed relative Intensities of the two data sets. An 
empiricle p-nrameter (♦) was defined to express the comparat ili ;;, of the 
.results, such that:- 

I' I" 2 

i ^ ) (6) 

uax *max 

^^ere n is the pointing location identifier. I' and I" are the intensities 

recorded by the rocket and ATM instruoents in their own respective units. 

I' and I" are the values of the intensities of the brightest point in 
max max 

the pattern. We found that there was a neglibable translational motion and 
a small rotation ('V' 2^) required to find the minimum of ^ within the uncer- 
tainty of the rocket pointing locations. The rotation is of little consequence 
to McMath 12379 as it is near to the disk center but could have been important 
to McMath 12387 which was on the eastern limb of the sun. The coordinates 
of the sun center and the active regions could be derived by overlapping a 
Stoney-Hurst disk, Figure 6. 

(iii) Analysis of the SO-82A Data. 

Two exposures were available of these regions (1A233 - a short 
wavelength plat<2 -9.5s; and l/v235 - a long wavelength plate -79.5s). A total 
of 10 scans of the entire active regions were made of the Images at various 
wavelengths using the microdensitometer at NRL. Due to calibration problems 
some of the low temperature lines of oxygen were not available to us. However, 
scans were made of both regions in Ne VII (465X), Si XII (A99%); Mg IX (368%), 
Fe XV (284^) and Fe XVI (33SX or 361%). 

Solving the coalignment problem was achieved by measuring the 
position of the disk center of the He II (304^) and Ne VII images from the 
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coordinates (in the mlcrodensitooeter's system) of 4 orthogonal limbs. 

Knowing the dispersion of the plates the position of the disk center could be 
derived for every other image. The orientation of the image is also known, 
hence the positions of the rocket experiment's fields of view can be cal- 
culated. See Figures 7(a) and 7(b). 

We derived the equivalent EUV line intensity for each of the above 
lines at every stable pointing position by convolving the intensity of each 
resolution element (1 arc sec) with the collimator efficiency at that position 
and summing the results. The SO-82A calibration curves were supplied by 
K. Dere. 

(iv) Combining the X-ray and ELT data to produce differential Emission 
Measure Models. 

Again we used Withbroe's (1975) technique to model the differential 
emission measure as a function of temperature. The data is presented in Table 
V. The quantity in parenthesis is the "fitting ratio" and is defined as 
the ratio of the observed to predicted (from the model) line intensities. 
Generally, the data seemed internally consistent with the lines conforming 
to the models well. Three lines were chosen in each model to have the heaviest 
weight in the fitting procedure. These were the brightest lines that defined 
the full temperature range (Ne VII, 0 VII, Fe XVII). The models are shown in 
Figures 8(a) and (b). The integral emission measure is also listed on the 
figures. The models were stable, with the exceptions of McMath 12379-7, 

Move to 12387, McMath 12387-10 and 13, which are not shown in the figures. 

The models have several interesting characteristics. The low 
ten^erature end is defined by Ne VII which gives rise to a plateau, generally 
there is a small peak in the distribution at about 2-3.10^ followed by a 
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sharp downward slope at higher temperatures. However, some distributions have 
no peak and the emission measure drops away from the Ne VII defined plateau. 

The significance of the gradients and maxima are reduced by the absolute 
intensity uncertainties of the lines ( 40%) and the pointing errors. However, 

the emission measure peaks bear a close similarity to those published 
by Withbroe (1975) and Levine and Pye (1979). 

The major anomaly in the data is the large discrepancy between the 
models and the intensity of Fe XV and Fe XVI. Inclusion of these data in 
the analysis never led to a stable model, hence they were excluded. The 
source of this anomaly has not been determined as yet but could be due to 
blending in these images with other features. 

(v) Comparison of the Emission Measures derived from different sources. 

SO-54 also had good coverage of McMath 12379. These data were 
analyzed to deduce the intensity and emission measure contours for this 
region. See Figures 9 and 10. Using the pointing positions established from 
the ATM orientation and measurements of the white light disk center, we 
derived the approximate emission measure from these contours. The results 
are listed in Table 6. 

The results indicate that the emission, as seen by the X-ray 
spectrometer and EUV spectroheliograph, are consistent within measurement 
uncertainties where as the SO-54 data consistently gave a low value in the 
emission measure by a factor of between 2 and 7. The source of this dis- 
crepancy is as yet undermined. We estimate that the total emission measure 
for McMath 12379 is 5.10^^, compared to 8.10^^ for McMath 12387. 
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TAB!,E V 

Modelling Parameters (ergs/cm^/sec/st) 
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Not used In modelling procedure as fluxes anomlnously high. 
Fit forced on these lines. 
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TABLE III 


Pointing Position 

Observed Emiss 

A 

ion Measures 


X-Ray Spectrometer 

ELW X-Ray ‘ 

S054 

McMath 12379-1 

1(47) 

2(47) 

5(46) 

McMath 12379-2 

2(47) 

2(47) 

8(46) 

McMath 12379- 364 

2(47) 

2(47) 

3(46) 

McMath 12379-5 

5(46) 

4(46) 

6(45) 

McMath 12379-6 

5(46) 

4(46) 

6(45) 

McMath 12379-7 

7(46) 

No Fit 

8(45) 

McMath 12379-8 

7(46) 

8(46) j 

3(46) 


* 


1.10 
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(d) Analysis of McMach 12664 


The philosophy behind this analysis is basically similar to that applied 
to the 11th June 1973 data. Table IV lists the available data. Again 
resolving the aspect was the major problem in the analysis. 

(i) Aspect Problem. 

We used two methods to try to establish the relationship between 
the rocket's aspect solutions and the images from SO-54 and S055. 

The pointing positions' coordinates of Aerobee 26.003 were derived 
as for 26.002. The problem reduced to relating the X-ray filtergram and EUV 
rasters to these positions. 

The first method (Figure 11) we tried utilized the position of the 
major sunspot in McMath 12664 as a reference point. The SO-54 filtergram 
had a corresponding white light image that could be pin registered, so that 
the X-ray data can be directly aligned to the svinspot. SO-55 then could 
be registered to this using EU\’ images taken in lines where this feature 
could be identified (C III) and its known rotational offset. 

The second method (Figure 12) involved using the coordinates of 
the two sunspots as measured by R. Levine from a magnetogram and marking 
these to scale with the rocket pointing posit.lons on drawing of the solar 
disk. A fit between the limb and the sunspots were then performed using the 
SO-54 white light image. Thus, we established the coalignment of the X-ray 
image and the spectrometer data. The SO-55 raster was aligned as before. 

This second method proved to be much more successful. The first 
method gave an unacceptable discrepancy between the white light limb and the 
drawing (solid and broken lines) as well as a poor correspondence between 
the X-ray line intensities as observed by 26.003 and the position of the 
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x-ray emission in the SO-54 filtergrams. We estimate the additional un- 
certainty introduced by this procedure to be about + 20 arc seconds. 

(ii) EUV and X-ray line intensities. 

The collimator functions were then convolved with the EUV pixel 
intensities at the resulting pointing positions, thus deriving the equivalent 
EUV line intensities. 

Table 6 lists the X-ray line intensities from Aerobee 26.003 which 
observed McMath 12664. These were to be modelled at HCO by Levine, however, he 
has found an anomaly between the data sets. Apparently there is a factor of 
3-5 excess flux in the X-ray line intensities over the level to be compatible 
with their EUV lines measurements. This problem could be due to a cross 
calibration error or a mismatch of the pointing positions. This is not 
resolvable but the most likely explanation lies in the 30 arc second un- 
certainty in the relative pointing solutions. Modelling such mismatched data 
would not be useful. It is interesting to note that the X-ray line intensi- 
ties from 26.002 seemed low with respect to NRL*s data but seem high for the 
second flight of the same instrument. This implies that either that there 
is a significant change in the X-ray flux for a relatively small change in 
the EUV flux or there is a cross calibration problem between SO-82A and SO-55. 
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GRAPH 1: Sunspot Alignment Procedure (^iethod I) 
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23 January l979 


GRAPH 2: Aerobee 26.003 Aspect Solutions (Method 2) 
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Conclu d ing Remarks 


During the course of this research project we have attempted to codify 
several different data sets to produce models of the emission measure dis- 
tribution in three active regions. In conducting this investigation we have 
had to develope methods of equivalencing data from X-ray spectrometers. X-ray 
filtergrams and EUV images. The primary advantage of adopting this technique 
is that it extends the temperature coverage of a model beyond that achievable 
by any single source of data. Where the temperature regimes are duplicated 
the secondary data may be used to confirm the model, as with the filtergram 
data (5054 and 5056). If there were an anomally then these data could 
indicate the presence of higher temperature material or a cross calibration 
error between the experiments. 

The models that were produced using 50-82A data with the X-ray line 
intensities from the sounding rockets demonstrated the feasability of this 
technique. However, we identified some areas that require further consideration 
if this type of work is to be attempted in the future, as will inevitably 
be the case with 5MM. They aret- 
(i) Instrumental Calibration. 

In order to put weight on any intensity measurements from a given 
experiment, a thorough knowledge of the instrumental sensitivity must be 
established prior to launch. Uncertainties of 40% generally render the results 
of differential emission measure modelling to be unusable for accurate physical 
interpretation. It is possible to establish the interrelation of the ATM 
instrument calibrations with a complete set of observations on a single stable 
region by all four instruments (5054, 5055, 5056, 5082A) . This has not been 
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done to our knowledge. Unfortunately the data available at the tine of our 
flights is not complete otherwise this would have been an ideal opportunity 
to attempt this. 

(ii) Coalignment. 

As the spatial resolution of Bragg Crystal Spectrometers improves, 
the ability to reconstruct the alignment of one instrument to another becomes 
a more important factor in a successful analysis. Hence, a reliable system 
to reproduce the relative position of a given resolution element (pixel) in 
one image to that in any other is vital. Tills was not available to use in 
this investigation. With our field of view ('v 35" FWHM) we were able to begin 
to derive the Intensity variations within the regions that were observed. This 
required accurate knowledge of the relation between the rocket pointing posi- 
tions and the aspect information for each experiment. 

(iii) Atomic Parameters. 

This investigation also hilighted the need for a uniform set of 
atomic parameters such as ionization equilibrium. Often we had to use several 
different sources of the ionization equilibrium to accomodate all the lines 
used in the analyses. We generally used Jordan’s calculations (1979) but 
used Davis et al (1977) for Fe. Elemental abundances also are a problem area. 

In our X-ray data we found the relative abundance of 0:Ne to be about 6.0 which 
is significantly lower than the value adopted usually (Withbroe, 1975). As these 
observations agree with other recent X-ray investigations (Parkinson 1978; 

Davis et al., 1975; Acton et al., 1976), yet are very different from the E^JV 
results, this may indicate a real variation in the abundances of one or both 
of these elements in these two temperature regions of the solar atmosphere or 
a problem in the Interpretation of the data. These aspects of the abundance 
problem merit a more detailed study, such as will be afforded by SMM. 
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(iv) Spatially Resolved Data. 


It was very useful to have available to us in the data analysis and 
interpretation phases simultaneous X-ray and EUV imaging 
information. This is particularly the case in the analysis of the resonance 
scattering problem where the geometry of the region was important. 

As a result of this man-year of research effort, we have produced 
draft copies of two papers to be submitted, as well as a proposal to study 
the resonance scattering problem and its implications in detail. We have 
developed programs to implement the analysis and interpretation of several 
different data types, which are summarized in Appendix IV. 


Acknowledgements 


The authors would like to thank all those that aided us on this Investi- 
gation. Our work with the SO-54 data at the Marshall Space Flight Center In 
Huntsville, Alabama was assisted by Ed Reichmann and Bill Henze, who helped 
us correlate the data and did the analysis on IDAPS. At the Naval Research 
Laboratory In Washington D.C. we were helped by G. Brueckner and Ken Dere. 

At Harvard College Observatory the data was processed by Randy Levine and 
George Withbroe who also advised on the differential emission measure modelling 
techniques. We worked with M. Gerrasimenko on the SO-54 data for McMath 12379, 
who processed the images and did the reduction of the data to derive Intensity, 
temperature and emission measure contours for this region. Dave Webb and 
A1 Krieger worked on McMath 12664 at American Science and Engineering with us. 

We would also like to thank the Lockheed personnel who helped with this 
project. John Kulander was involved closely with the development and execu- 
tion of the res^nance scattering study, and responsible for developing the 
theory. Tim Roethig has assisted us with advice on computing and made the 
HPIOOO computer system available to us. Various men±>ers of the scientific 
staff at LPARL have been given useful comments during the course of this 
project . 


44 


Figure Captions 


Figure 1: 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Six typical spectra from the collimated Bragg crystal spectro- 
meters flown on Aerobee 26.002 and 26.003. 

Two typical fits of the Lorentzlan profiles compared to our data 
for Fe XVII scans. The raw data are designated by a series of 
points and the fitted spectra by a solid curve. The numbers 
indicate the positions of stable fitted lines in each spectrum. 
In each case 30 prospective lines were chosen arbitrarily to 
represent the spectra. 

Analysis of the 0:Ne abundance ratio. The histogram represents 
the frequency of occurrence of a given line ratio. The family 
of curves represent the variation of this ratio as a function of 
electron temperature for various abundance ratios of 0:Ne. 0:Ne 
ratio * 6(^.5) most closely reproduces the first peak in the 
distribution at the measured coronal temperatures. The second 
peak in the distribution can be explained in terms of a higher 
relative abundance ratio or a considerably higher temperature. 
Representative 0 VII spectra from each of the four types of 
regions observed by the two rockets. The Forbidden and Inter- 
system lines are all drawn the same size thus demonstrating the 
change in the resonance line flux. 

X-ray line intensity maps for the Regions: McMath 12379, 12387 
and 12664. The lower temperature lines show less contrast but 
more extensive emission than the higher temperature lines. They 
are compared to two intensity maps from SO-56 drawn to the same 
scale. 
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Figure 6(a)&(b): Rocket pointing locations with respect to solar disk for 

both flights. 

Figure 6(c): Shows a Stoneyhurst overlay of an SO-54 image used in co- 

alignment procedure for Aerobee 26.002 (courtesy of MSFC) . 

Figure 7: Sequence of photographs of (a) McMath 12379 and (b) McMath 12387 

in different temperature lines compared to an X-ray filtergram 
and the pointing locations of the rocket fields of view (courtesy 
of NRL and MSFC). 

Figure 8(a)&(b) Differential emission measure models derived for each 

field of view where a stable model was found (a) McMath 12379, 

(b) McMath 12387. The models are calculated by combining the 
X-ray emission line data from the rocket flight and the EU\’ image 
data from SO-82. 

Figure 9: Intensity map of McMath 12379 derived from SO-54 filtergram 

(courtesy of AS&E) . 

Figure 10: Emission measure map of McMath 12387 derived froma ratioing 

of two filtergrams from SO-54 (courtesy of AS&E). 

Figures 11 & 12: Relative pointing positions of SO-54, SO-55 and the rocket 

fields of view as derived by two different coalignment techniques. 
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A ctive Regions from Observations of Coronal X-Ray Emission Lines 


by 
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Abstract 


Observations of three solar active regions, obtained by rocket borne soft 
X-ray Bragg crystal spectrometers, are presented. The analytical methods used 
in the analysis of these data are discussed and the results presented. We 
determine various parameters which describe the state of coronal plasma at the 
time of the observations in these active regions. The anomalous behavior 
of certain line ratios of helium-like ions give strong evidence for resonance 
scattering (f lourescence) of X-radiation in the corona. During this work 
we have derived the relative abundances of some elements, we describe in some 
detail the results for the 0:Ne abundance ratio. 
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I. Introduction 


The soft X-ray spectral band (l-3oX) of the solar spectrum has a pre- 
dominance of emission lines that mainly originate from coronal Ions in the 

hydrogen, helium and neon isoel ectronic sequences. These emission lines are 

8 10 “3 

characteristic of a low density plasma (10 - 10 cm ) maintained at a high 

temperature (10^ - 10^ ^K) • The line intensities originating from these ions 
contain information on the prevailing physical conditions of the plasma and 
quantity of emitting material in the observed region. 

The work presented here describes a series of observations of various 
coronal features obtained during two flights of an array of six Bragg Crystal 
Spectrometers. These instruments monitored most of the prominent emission lines 
betw'een 9S and 22. sS. Data were gathered on three active regions and a region 
of quiet corona. The field of view of the spectrometers were restricted 
in order to isolate the emission from, a descrete region. The spatial resolu- 
tion varied from about 35 arc secs to 150 arc sec. full width half maximum. 
(FVHl'O . The experiment package included an interactive pointing control 
system which enabled the experimenters to perform, raster scans over any regions 
of interest . 

The following sections contain a description of the payloads, the data and 
the techniques used in the analysis. We have derived the line intensities of 
all the X-ray emission lines observed in the spectral ranges covered by the 
spectrometers for each stable pointing position. For reasons of brevity ve will 
only give the mean value of any derived quantity for each distinct coronal 
region where tabulating the full set of results would not add to the informa- 
tion presented by this work. However, individual results that are especially 
interesting or anomalous will be discussed rts they arise. Where possible, we 
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have estimated the uncertainties on any derived quantity. We have assumed 
pesamistlc values for effects of possible s'^urces of error during these cal- 
culations. Some indication of the conservative values put on .these uncertainties 
is given when we have several measurements of the same quantity, generally the 
values fall well within the calculated- uncertainty . In Section IV we discuss 
the results of the analysis and their interpretation. In the last section 
the limitations of our data and its possible uses are discussed. 

II. Instrumentation and Observations 

The experiment consisted of an array of six scanning Bragg crystal 
spectrometers and a bent crystal spectrometer (Acton et al.,1976). Each 
spectrometer monitored a limited wavelength range centered on a selected 
prominent emission line in the soft X-ray spectrum. The field of view of each 
spectrometer system was restricted by a two dimensional multigrid mechanical 
collimator. The ability to limit the field of view, coupled with a SPARCS 
inflight command link, made it possible to execute raste- maps of regions of 
interest under real-time control of the experimenters. 

Table 1 lists the primary characteristics of the payloads for the two 
flights. On the first flight (Aerobee 26.002) of this payload four of the 
Bragg Crystal spectrometers had 35 arc second (FWHM) fields of view whereas 
the other two instruments had 150 arc second fields of view. These instruments 
with their poorer spatial resolution duplicated two of the higher resolution 
instruments, and were designed to monitor the active region as a whole. They 
had better counting statistics, owing to the Increased signal levels, 
than the 35 arc second systems. They were replaced on the second flight 
(Aerobee 25.003) by instruments that monitored other spectral regions at the 
higher spatial resolution. 


L 
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The collimators' fields of view were circular and had a radial sensi- 
tivity dependence represented by the function, P(r), given by: 

P(r) « 2(<f - r sin (f)A, (1) 

where r is the fraction of the total radius of the field of view (r *= 1 f 

max 

5/A FWHM) . ({) is the angle who’s cosine is r. The collimator transmissions were 

measured before each flight. The offset of the first side lobes were set to be greater 
than 32 arc mins from the optical axis to prevent contamination of the data from 
other active region emission. 

The diffraction cr^'stals were KAP for all the spectromters, except the Mg XI 
system flcwn on Aerobee 25.003, which used ADP. The crystals were mounted on 
platforms that were scanned smoothly and continuously through an angular range 
of about 3^ on a motor-driven cam. The crsytals were offset with respect to 
one another to access the appropriate angular (wavelength) ranges. They were 
grouped in two sets of three,, each set with an independent drive system and 
shaft encoder. The encoder divided each scan cycle into 512 angular bins. 

A scan over the total spectral range was completed for both the increasing and 
decreasing wavelength directions in about 14 seconds. The integrated reflec- 
tivity of the crystals were determined from measurements made by R.L. Blake 
and A.R. Burek. O 

The proportional counter detectors were filled with an argon-methane 
mixture at a pressure of 1 atm. The polypropylene windows were carbon coated 
and were 1 pm thick. The number of events registered by the proportional 
counters were accumulated in registers before being transmitted to the ground 
system on the P.C.K. telemetry. The registers were scanned at a rate of 256 Hz. 
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The SPARCS inflight-conunand link was used by the experimenters to set and 
maintain the rockets' pointing with respect to the solar regions of interest. 

A diffraction pattern, produced by special grids in the collimator stack, was 
superimposed on a hydrogen-alpha solar image produced by a telescope/filter 
system. One of the spots of the pattern was designed to correspond to the 
optical axis of the collimation system. The image was fed into a television 
camera and transmitted to the control center to aid the real-time control of 
the pointing position. The image was also recorded by a 16 mm film camera in 
the experiment package for post flight analysis. 

A brief description of the experiment has also beein given by Acton and 
Catura (1976). Table II details the observations made by the instruments. Most 
of the observations were coordinated with the ATM experiments. Figures 1(a) - 
(f) show typical examples cf the spectra obtained by these experiments. 

III. Analysis 

At every stable pointing position the data accumulated by each of the 
six spectrometers consisted of a pair of spectra, scanned in the forward 
(decreasing wavelength) and reverse directions. A spectrum consisted the 
number of events recorded at each crystal encoder position. Due to non 
linearities in the crystal drive cams the count rates had to be normalized over 
the observation time of that position. In order to interpret these spectra, the X-ray 
lines recorded in the data have to be identified and their intensities derived for 
each field of view. These results then have to be coordinated with the aspect 
information and any data gathered simultaneously on the same regions by 
other instruments. 
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The spectral regions, observed by these instruments, have been extensively 
monitored by previous experiments (Parkinson, 1975; 

Acton et al., 1976). Hence, the wavelength of the prominent lines are 
familiar to us. The data were of sufficient quality to be able to resolve 
many of these lines. All the spectra were analyzed by fitting the following 
expression for the number of events (K) in an encoder bin (i) over the complete 
spectrum. 


N. 

1 


n=L 

B + I. 
n=l 


An 


1 -f A 


cts/sec/encoder bins. 


( 2 ) 


where B is the background count rate in each data bin, which is assumed to 
be constant over the whole spectrum. L is the total nuniber of lines present 
in the spectrum. is the peak amplitude of the nth line and is given by: 
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S 

n 


a - X.) 

n 1 


(3) 


where w is the full width at half maximum of the nth line and X is its 
n n 

position. is the position of the ith encoder bin. X^ and 

are expressed in the same units (A, i) . The line profile ex^pressed in Equation 
2 takes the form of a Lorentzian, which we found to be an acceptable approxi- 
mation to the observed line profile, which in reality is a convolution of the 
natural doppler profile, the collimator function (Equation 1) and the crystal 
rocking curve. The amplitude, position and half v;idth of each line, as w^ell 
as the background, were free parameters in the least squares fitting calcula- 
tion. There were consequently a maximum of 3n 1 parameters in the fit. 
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As mentioned above, the positi' s of the brighter lines are well 
established in these spectral regions, also the characteristics of the flight 
crystals were well knonw. Hence, the derived positions of the lines and their 
half widths could be compared with these established results. The agreement 
was, generally, very good, giving us some confidence as to the integrity of 
the analysis procedure. Where fainter lines were blended but their positions 
are known from other observation.^ made at higher spectral resolution, we could 
just fit for the line amplitudes using a half width derived from the surrounding 
lines for the blends. 

The total number of counts observed in each X-ray emission line in unit 
time was derived by integrating over the fitted profile. This can be done 
analytically, resulting in:- 


c 

n 


ttA w 
n n 


(4) 


where c^ is the integrated count rate. We averaged the results from the 
forward and reverse scans of a given spectral scan for each pointing position 
at this stage. Where there was a large descrepancy in peak counting rate, 
line width or wavelength between the forward and reverse scans the data were 
considered to be invalid. These are elliminated lines that were miss matched 
or merely noise spikes. The uncertainties derived from the standard devia- 
tions on the fitted parameters were added quadratically , The total photon 
flux (i) could then be derived from the efficiencies of the instruments 
components : 

C a' 

* * "a "i” — photons/cm‘’/sec (5) 

eff c n 
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where u. is the average crystal scai; i.ue (rads/sec) and is the integrated 

reflectivity of the crystal (rads/sec), n is the collimator transmission, and e 

is the quantum efficiency of the proportional counters. The effective area of 

the instrument (A is derived from the geometrical characteristics of the 
et r 

individual instruments. These were determined by measurement and expressed 
by the empiricle relation: 

where the C’s are constants, deduced from the preflight calibrations, 0 

n 

is the Bragg angle of the line and is the wavelength of the line. 

Initially, the fitting was applied to spectra tl^at were summed over each 
active region. Thus, we were able to establish the mean positions, widths 
and intensity of the lines to use as first approximations for fits to the 
brighter lines in the individual spectra. The intensity of all lines identi- 
fiable in the spectra from each stable pointing position were analyzed as 
described above. This technique makes it possible to analyze a spectrum as 
a whole, taking into account contributions to the background of the wings of 
all lines in the spectrum. Failure to do this will lead to an over- estimation 
of the background signal, and consequently an underestimation of the line 
intensity. For example, by fitting a triangular profile to a line giving 
it the same FWHM can lead to a possible 30rc under estimation of the integrated 
line intensity, from neglecting the wings of the line, assuming the background 
is perfectly determined. 
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IV. Results 


(a) Emission Line Intensities 

There were a total of 31 stable pointing locations from the two 
flights. Each location yielded a total of 12 spectra. Using the above 
techniques we positively identified 58 emission features in the 9 - 23% 
band of these spectra. Table III lists the wavelength and possible identi- 
fication of these features. The wavelengths derived from each flight ate 
listed separately .with the estimated uncertainty in parenthises, and a 
weight mean from these values calculated. These are compared to the results 
from other observational and theoretical investigations. Generally the 
agreement is very good. 

The volume of results from the spectral fitting of the data from 
the two flights preclude their inclusion in this work. However, the mean 
intensities of the lines are listed in Table IV for each of the three active 
regions and the quiet coronal region for the high and low spatial resolution 
data. The values in parentheses are the factors between the brightest and 
faintest lines recorded for that region. Upper limits on the intensities are 
put on lines where no fit was obtained. 

These data show the contrast between the quiet corona and active 
region spectra. In the ’’cooler** lines, such as 0 VII, the quiet corona 
registers a significant signal, although it is a factor of 6 to 10 weaker 
than the corresponding active region spectra. However, all other lines formed 
at higher temperatures show little evidence of emission from the quiet corona, 
within the sensitivity range of this experiment. Also the contrast between 
the various active regions is evident from these data. In all lines, McMath 
12664 is the brightest region by a factor of about 3 to 5 over McMath 12379 
and McMath 12387 which were of comparable brightness. 


The brightness factor results demonstrate the advantages of having 


higher spatial resolution available in such an experimental package. The 
130' arc second FVHK field of view instruments all gave brightness factors 
of between 1 to 4, indicating that the instruments were incapable of giving 
any useful information about the structure within an active region. However, 
the higher resolution data gave brightness factors of between 3 and 30. 

Thus it is possible to use these data to produce intensity maps of X-ray 
regions observed. The results of the mapping procedure can be seen in Figure 2 
for the three active regions observed. Each stable pointing position is 
assigned a square with a side of twice FWHM. These are shaded in accordance 
with the relative brightness of the primary lines observed by the flights 
(0 VII - 21.6; 0 VIII - 18.95; Fe XVII - 13.01; Ne IX - 13.44; Ne X - 12.13; 

Mg XI - 9.17). The position with brightest intensity of any given linewas assigned lOi 
all other brightnesses were blocked into 3 groups in relation to the brightest 
position. The resulting maps are compared to an intensity contour plot from 
the S054 instrument on ATM (Gerassimanko, 1978) for two of the regions. The 
results correspond very well. One interesting feature of these plots is that 
the contrast between the "cooler" lines (0 VII and 0 VIII) is less than those 
of the other lines, the emission from which seems to be more localized. 


(b) Resonance Line Ratios 

The intensity of an X-ray emission line can be expressed in the 
form (Pottash, 1964) : 


E 


■D g f A X 



2 

G(T) N 

e 


dV, 
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where E is the line intensity (photo, .'cm /sec/st) at the Earth. D is a 
constant, g is the Gaunt factor and f is the oscillator strength. A is the 
elemental abundance and A is the wavelength (^) . is the electron density 
and the temperature function, G(T), is given by: 

. 5040 U\ 

1/2 T ' 

G(T) ■ T ' 10 N^, (8) 

where is the number density of the ion stage responsible for the transition 

at the temperature, T. W is the excitation energy of the transition. 

For any given transition equation (7) may be rewritten, assuming 

that the t-mission can be characterised by an isothermal plasma and 

2 

a suitable emission measure N V, 

e ’ 

E^ = D f^ g A^ Gj^ ^’^eff^ photon/cm^/sec/st (9) 


Although using this expression directly with a measured line intensity 
can yield information on the coronal plasma, such as its effective temperature 
or emission measure, the results are subject to uncertainties in the data. 

By ratioing line intensities such uncertainties can be minimized if care 
is taken in the choice of the lines. Hense the ratio of any two lines can 
be given by: 






5040 (W-W) 

-( i — — 

T 

10 eff 


( 10 ) 


It is evident that by choosing lines carefully some of these factors will 
be eliminated. Further, if the lines are observed by the same instrument 
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then the ratio will be free from calibration and instrumental uncertainties. 


For example, choosing lines produced by the same ion will give a ratio 
that is mainly dependent on the relative oscillator strengths. Table I 
lists oscillator strengths derived in this fashion. Similarly, it is possible 
to derive reliable relative abundances by choosing lines with similar 
temperature functions, and assuming their oscillator strengths (Acton and 
Catura, 1975; Patkinson, 1977). We investigated the relative abundances of 
Mg to Fe and Ne to l!g in this manner. 

The Mg:Fe ratio utilized the Mg XI line at 9.17A and Fe XV’II line 
at I 5 .OO 9 X. The relative abundance ratio was found to be 0.71(+ 0.15) from 
the observations of McMath 12664. This compares with values of 0.6 (Potash, 
1967), 0.97 (Parkinson, 1977) and 1.2 (Withbroe, 1971). The relative 
abundance of Ne to Mg was determined using the Mg XI line at 9.17^ and the 
Ne X line at 12.134.?. It was found to be 1.6 (+ 0.8) also from the observa- 
tions of McMath 12664. This value compares with 1.2 (Pottash, 1967), and 
0.07 (Parkinson, 1977) and 0.85 (Withbroe, 1971). 

The relative abundance of 0 to Ne has been the subject of several 
previous investigations. Its value seems to differ depending on the origin 
of the data. EUV data tend to give a high value for the 0:Ne ratio, for 
example, Withbroe (1971) gives the ratio to be 15.8. However, soft X-ray 
emission line data generally give lower values. Acton et al. (1975) have 
calculated a value of 4.7, Parkinson (1977) gave a value of 6.6. We use the 
same technique as Acton et al. (1975), which utilizes the intensity ratio of 
the 0 VIII line (18.95?) to the Ne IX line (13.445). Figure 3 shows the 
frequency distribution of the 0 VIIIrNe IX line intensity ratio from the three 
active regions observed by the rocket experiments. The ratio was calculated 
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with its uncertainty for each stable pointing position. Each value was then 
considered to be a nonnal distribution with a width equal to the calculated 
uncertainty and an integral value of unity. Tnese were then summed to produce 
the frequency distribution. This ensured that any ratios with a large 
uncertainty had a minimal effect on the results. Plotted above the resulting 
histogram are a series of curves that represent the temperature dependence 
of the 0:Ne ratio calculated for several values of this ratio. 

The observed frequency distribution has a distinctly dual population. 
The first with a maximum at a ratio of about 7.5 and the second at about 15. 

The second peak is mainly due to results from the brightest fields of view 
of Mcllath 12664. 

The effective temperatures of the regions were determined using 
line ratios of adjacent ionization stages of several elements. The results 
are shown in Table V. The particular 0 V1I:0 VIII ratio used xv’as chosen as 
the data were derived from the same instrument, as was the Fe X\^II:Fe XVIII 
ratio. Although the 0 VII resonance line at 21.602A is brighter and vrould 
give a statistically more reliable result, there is some evidence for a 
change of the resonance line intensity due to scattering so it was not used 
in this calculation. The evidence for resonance scattering discussed in a 
subsequent section. The effective temperatures determined from the line 
ratioing technique all lie in the 2.9 - 3.9.10^ range. Using these 
temperatures and the main peak of the frequency distribution, we 
derive an 0:Ne abundance ratio of 6. 0(^.5). This is consistent with 
Parkinson’s determination (1977). 

The second peak in the ratio frequency distribution would indicate 
a relative abundance of 12 (+2), which is more in agreement with the EUV 
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results, assuming the same temperature ranges. However, it is possible that 

we observed some of the flare associated plasma, hence an abundance ratio 

6 o 

of 6.0 would be obtained by a temperature of approximately 9.10 K. Non 
uniform abunadnaces throughout the solar atmosphere would be a difficult 
problem to explain and model. A further point of interest in this regard 
arises from the McMath 12387 data, where the brightest field of view gave a 
value of 2.5 for the 0:Ne line ratio, being responsible for most of the power 
in the frequency distribution below 4. This result cannot be lightly dis- 
counted, as it is not explainable with an abundance ratio of greater than 4 
at any temperature, and statistically significant. 

(c) He-like Ion Line Ratios 

There are two line ratios that are of particular interest in the 
data from the He-like ionization stages. Adopting the nomenclature of 

Gabriel and Jordon (1969) , R is defined as the ratio of the Forbidden line 

2 13 2 13 

(Is S - ls2s S) to the Intersystem line (Is S - ls2p P) . This ratio is 

9 -3 

sensitive to electron density above a threshold of 6.10 cm (Freeman et al., 

3 3 

1971) due to the possibility of a transition between the P and S levels. 

The possible temperature dependence of R has been investigated by several 

groups (Acton et al., 1972; Blumenthal et al., 1972; Gabriel and Jordon, 

1973). We are unable to determine the nature of any temperature dependence 

from these data within the experimental uncertainties. Hence, the ratio is not 

strongly temperature dependent which favours the views of Gabriel and Jordon 

(1973). The density for all our regions lies near or below the threshold 
9 -3 

value of 6.10 cm 
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The second ratio is designated as G and defined as the sum of the 


intensities of the Forbidden and Intersystem lines divided by the intensity 

2 1 1 

of the resonance line (Is S - ls2p P) . This ratio is independent of 
electron density as it contains the sum of the forbidden triplet transitions. The 
temperature dependence of the ratio has been studied by Acton and Brown (197S), 
and Meve et al. (1977a, b). 

The results for all the observed Helium-like ions are listed in 
Table VI for each region and spatial resolution. The G ratio show’s an 
interesting variation over different regions for 0 VII. These variations are 
significant with respect to the experimental uncertainties. The quiet corona 
gives the lowest value of G. From the theoretical calculations of Acton et al 
(1978) and Mewe et al. (1977) this trend is contrary to that anticipated from 
a low temperature region (£ 2.10^ °K) such as the quiet corona. Further, 
the active region results show significant variations that cannot be due to 
temperature differences alone. There is no evidence in the data for a 
significant departure from equalibrium conditions in the plasma that could 
explain the discrepancies between the calculated value of G (^/ 0.8) and the 
observed values. There is also an indication of a systematic difference 
between the results from the two different spatial resolutions from any 
given region. None of the other He-like ions show similar variations. 

We believe these effects are caused by resonance scattering of 0 VII 
in the corona. Recent work by Acton (1978) show that resonance scattering 
is a viable mechanism in the corona. The most susceptable ion is to this 
effect at normal coronal temperatures (^ 3.10^) is 0 VII. The variations 
observed in the G ratio listed in this work could be accounted for by an 
enhancement or a reduction in the resonance line intensity due to scattering 
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(Strong, 1978; Strong and Kulander, 1979). The forbidden triplet transition 
will remain unaffected by scattering effects. Hence, G will be Increased 
by a net absorption of resonance line radiation along the llne.of sight, 
such as was observed in the three active regions, or will be decreased by 
an enhancement due to scattering into the line of sight as observed in the 
quiet coronal observations. 

(d) Modelling of the Active Region Plasma 

The X-ray line intensities may be used to derive differential 
emission measure models of the plasma in the observed regions. The resonance 
lines of fine ions were used in the modelling from the 35 arc sec FWHM fields 
of view, the parameters are listed in Table VII. The modelling of the dif- 
ferential emission measure involves the use of absolute line intensities. 

These are subject to statistical fluctuations in the observed count rates, 
which can be quantified, and also to errors in the relative instrumental 
calibrations. Any atomic parameters used in such a calculation are also 
subject to errors. Every effort was made to minimize the relative preflight 
calibration errors in the spectrometers, however, some residual differences 
will remain when comparing the fluxes observed by different detection systems. 
The most important factor is the uncertainties introduced by the atomic data, 
such as the elemental abundances, oscillator strengths and ionization balance. 
While we used the best atomic data available, there still remains the 
possibility of introducing errors into the calculation from this source. 

The modelling techniques used in this investigation was first 
described by Withbroe (1975). This method compares the intensity of a number 
of emission lines as calculated from a differential emission measure model 


with the observed intensities. The discrepancy between ther. is used to, 
suitably weighted by its emission function, to adjust the model. The 
procedure is repeated until the model has converged to give th.e observed line 
intensities. The ability of this method to resolve structure in the dif- 
ferential emission measure model is dependent on the vidthsof the temperature 
functions and their relative positions in the temperature domain. Figure 4 
shows the temperature functions of the five lines used in the investigation. 

The temperature functions of X-ray emission lines are generally quite broad, 
and are consequently, incapable of defining small or narrow features in 
the derived differential emission measure function. 

The speed of convergence is dependent on the number of lines used 
in the modelling procedure, the amount of overlap between their emission 
functions and the internal consistency of the data. Although many lines 
were available to us in the analysis, many of them would be redundant and 
slow the convergence. Convergence was achieved, on average, in about 50 
iterations. The convergence was stable, and could be taken to several hundred 
iterations without significant change. As mentioned above the 0 VII line, 
which was the most sensitive to lower temperature material, was corrected 
for resonance scattering effects. Failure to do this led to an unstable 
result. The 0 VIII line may also be subject to this effect, to a lesser 
extent, but was not corrected. 

Table VIII lists the results of the modelling. Generally, the 
models could be characterized by a single peaked function. The table gives 
the temperature of the peak of the modelled distribution and the i/e width. 
The integral of the differential emission measure is calculated. The iteration 
procedure was terminated when all the line intensities lay within their 
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calculated unrerff^iiritv and no significant (< 1%) change occurred in the model 
over successive iterations. The ratio of the predicted to observed intensity 
was also calculated (the fitting ratio). 

The models indicate a mean value for the temperature of the regions 
similar to those derived from the line ratio technique described above. 

The brightness of the active regions does not seem to depend heavily on the 
effective temperature but more on the total emission measure. Hence, it would 
seem that the density and extent of the active regions determines the X-ray 
emission rather than their temperature. The 0 VII predicted line intensity 
seems to be consistently low compared to 0 VIII line intensity. This could 
be due to resonance scattering effects in 0 VIII that were not taken into 
account of the calculations or a further discrepancy in the 0:Ne abundance ratio. 

Generally, the data fitted to within 10-40% of -.he absolute line 
intensity. This is quite good considering the inherent problems in assigning an 
absolute line intensity as mentioned above. 
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V . Con cluding Remarks 


In this presentation we have described the reduction and analysis of the 
data accumulated during two flights of an instrument consisting of an array of 
six Bragg crystal spectrometers. This payload gathered spectra in the soft 
X-ray wavelength band (9-22S) from three active regions and an area of quiet 
corona. We have described a technique to reduce these data to line intensities 
that give a reliable result. We then demonstrate the variety of information 
that can be derived from such results while indicating their limitations and 
uncertainties . 

We find that the spectra should be fitted as a single entity, and not 
treated as a series of individual line profiles. For our instruments, with 
their relatively small solid angle fields of view and RAP crystals, we determine 
that a Lorentzian profile adequately fits the data. By fitting all the lines 
in a single spectrum simultaneously using a representative profile with an 
itterative least squares technique, we find that we get r more accurate 
determination of the background level (comprised of instrumental noise, X-ray 
continuum, and cosmic ray events). Thus, avoiding overestimate of the back- 
ground due to contributions from the wings of surrounding line profiles, and a 
consequent under estimate of the integrated line intensity. Errors can arise 
of over 10 % by considering the background level to be determined by a portion 
of a spectrum *'suf f iciently far from the line** - < 7 FWHM*s from the line center. 
Often there are no parts of the X-ray spectrum sufficiently devoid of lines, 
resolved or unresolved. 

An important result is the variation in the G ratio of 0 VII. We are able 
to show that these variations are significant with respect to the uncertainties 
in the data. We can not interpret the observed ratios in terms of temperture 
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variations from one region to another as the trend is contrary to that pre- 
dicted by recent theoretical studies (Acton & Brown, 1978; Mewe and Shrivjer, 
1977a, b) . Nor is non-equalibrium conditions in these regions a likely 
explanation as it firstly requires that the quiet corona is in a continual 
ionising state over a large area to account for its low G ratio, and that at 
all 26 stable pointing locations in the three active regions there was a pre- 
dominance of recombining plasma in the field of view. This would be an un- 
likely occurence. Also, none of the other He-like ions show similar variations, 
within the limitations of experimental uncertainties. We find that the data is 
interpreted most convincingly in terms of resonance scattering (flourescence) 
in the coronal plasma. These data have been investigated in detail as a result 
of this work (Strong, 1978; Acton, 1978; Strong and Kulander, 1979). This 
phenomenon could be utilized to study the plasma in more detail than has been 

possible previously as scattering effects can give us information on the density, 

and distribution of the coronal material. It should also be possible to detect 

the polarization of X-radiation if these conclusions are correct. 

We have determined the wavelengths of the prominent lines observed by 
these instruments, and proposed possible identification for them.. In taking 
the line intensity ratio of several lines we have determined several 
relative oscilator strengths and abundances. One particularly interesin,; 
result of this work is the relative abundance of 0:Ne. This statistical study 
gives a relative abundance ratio of 6.0(+0.5) for active region temperatures. 
However, there are indications of both significantly higher and lower values of 
the ratio from different regions. The value of 12 (+2) obtained from the 
brighter regions of McMath 12664 can be Interpreted alternatively in terms of 
a higher temperature plasma being observed whilst maintaining the lower abundance 
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ratio. There is some support for this view as some flaring activity was 
observed in this region just prior to the observations. However, it is dif- 
ficult to explain the lower value obtained from the brightest field of view of 
McMath 12387. The possibility of an anomalous distribution of 0 and Ne can be 
resolved by a further study of this problem by the XRP instrument on SM>'. 

We derive differential emission measure models of the coronal plasma in 
these activ ; regions using a technique pioneered by Withbroe (1975). We find 
that, for the most part, the emission measure distributions as indicated by 
the analysis can be represented by a single isothermal peak. This is indica- 
tive of the poor temperature range and resolution of the lines available to us 
in these data. The X-ray data is capable of representing the emission from 
the hot active region plasma, and indicates a uniformity of temperature through- 
out a region. This would be an expected result due to the high conduction of 
the coronal plasma. However, these data can give us little information on the 
lower temperature material at the base of these active region loops. This 
indicates the need of coordinating this type of data with EUV observations of 
the lower corona and transition region. This will be the subject of a future 
study. 

The type of data derived from soft X-ray emission line spectra are subject 
to many types of uncertainty, it is important to understand these in order to 
be able to know how to interpret these data reliably. The derived line 
intensity is subject to statistical limitations on the observed count rates. 

Also the preflight calibration of all the instrumental efficiencies is ver>' 
important. An accurate knowledge of these makes it posssible to calculate 
absolute line intensities of between 25^ and 40% uncertainty with these data. 
Assigning identifications to the spectral emission features is only a problem 
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for faint or blended lines, as most ol the positions of the prominent lines 
are well known and it is possible to derive the wavelength of the lines to 
an accuracy of ^ 0. oosX in general. 

The interpretation of these data requires the use of parameters such as 
abundances, ionization equalibrium calculations and various atomic parameters, 
all of which can introduce uncertainties into the calculations. We show that by 
using the technique of line ratioing it is possible to minimize the effect 
of instrumental calibration errors, and uncertainties in the atomic parameters, 
while still being able to derive useful characteristics of the particular 
ionization stage or the plasma that produced the emission. Consequently, in 
spite of the above difficulties, it is possible to derive new and interesting 
results from this type of X-ray data. 
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Table VII 


Ion 

0 VII 

0 vni 

Fe XVII 

Ne IX 

Mg XI 

11 

1 

2 

4 

3 

. 5 

X(A) 

21.60 

18.95 

15.01 

13.44 

9.17 

f 

0.69 

0.42 

2.2 

0.72 

0.75 

A 

2.4.10“^ 

2.A.10"^ 

5.2.10"^ 

4.0.10”'^ 

4.5.10"^ 


2.0 

3.2 

4.4 

3.6 

6.3 

iT(l/e) l.>3.3 

Pointina Position 

2. 0-5.7 

2. 5-6. 8 2. 1-5.9 

2 

Data (Dhts/cni /sec) 

3.5-10.0 

MM 12379 -1 

24300 

38000 

34000 

5350 

- 

-2 

45000 

62700 

52700 

12200 

- 

-3 

45500 

65200 

32200 

10200 

- 

-4 

38300 

71000 

35200 

7050 


“5 

15900 

8100 

2800 

- 

• 

-6 

9500 

11500 

7800 

500 

- 

-7 

13600 

14200 

5300 

2800 

- 

-8 

23700 

18500 

- 

1100 


MW 12387 -9 

36700 

17700 

4100 

1650 

- 

-11 

47500 

36700 

16000 

2500 

- 

-1? 

46000 

54300 

38300 

6300 

- 

-14 

58600 

88600 

104000 

9530 

— 

MM 12664-1 

274000 


352000 

88300 

9950 

-2 

177000 

- 

347000 

61300 

11100 

-3 

82000 

- 

157000 

25000 

6500 

-4 

34200 

- 

35200 

7350 

1000 

-5 

46500 

- 

29600 

3350 

1250 

-6 

120000 


90100 

12b00 
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Table XI 

26.003 Oxygen 0 VII & 0 VIII Line Intensity (Photons/cm^sec) 
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Probably affected by atmospheric absorption. 
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Abstract 


We (femonstrate that the line intensity ratios of certain He- like ions can 
ejshibit significant variations from the predicted values, and that these observed 
variations are most probably caused by the scattering of X-radiation in the 
corona. We develop the equations of radiative transfer for a sinple coronal 
model, and investigate the significance of scattering effects as a function of 
the basic model parameters such as density, teinperature , incident radiation field 
and orientation of the region to the observer's line of si^t. The model is 
applied to several of the observed situations using the 0 VII ion as an example, 
and the results are compared with the data. The agreement is good for the case 
of the quiet corona vhere the model is directly ^plicable to the 
observations. A discrepancy of 27% in the resonance flux is adequately accounted 
for using the irodel. The agreemerit between the model and the active region 
data is not so convincing. We discuss the reasons for this. 


1 . Introduc±ion 


Ihe highly ioru: 2 ed cjoronal plasira emits a spectrum that is dcaninated by 
soft X-ray anission lines. The most promirient .Unes originate from transitions 
of the H-liJce, He-like and Ne-like isoelectronic sequences of the heavier 
elements. The fhysics of the production of these X-ray lines is generally 
well established as the ions have a siJiple electronic structure, and the lines 
have been widely studied both observational ly and theoretically (Gabriel and 
Jordan (1969); l^e and Shiijver (1978a, b) ; Bhalla et al (1975); Gabriel and 
Jordan (1973) ; Evens and Pounds (1968) ; Parkinson (1971) ; Pari ;’jison (1975) ) . 

The intensity of such lines and their ratios have been used to derive the 
fundamental characteristics of the coronal plasma, however iiiplicit in these 
analyses is the assunption that the X-rays, once emitted, escape from the 
corona without further interaction with the plasma. That is, the optical cfepth 
of the coronal material is negligibly small for any given soft X-ray line. 

The work presented here challenges this assunption by shcving that results 
derived from observational data indicate variations in certain He-like ion 
ratios that are attributable to resonance scattering effects in the ambient corona. 

Acton et al (1976) have presented cbservaticnal evidence that the assunption 
of an optically thin corona is rot valid. A subsequent theoretical stuc^’ 

(Acton, 1978) has shewn, for some lines, that column densities exist in the 
corona capable of producing observable changes in some X-ray line intensities 
due to absorption and scattering along the line of sight. 

He-like ions are particularly useful for investigating resonance scattering 

2 1 1 

as the resonance line (Is S - ls2p P) intensity can be compared directly 

2 13 2 1 

to the sum of the forbidden (Is S - ls2s S) and the intersystem (Is S - 

ls2p P) lines. The triplet lines are forbidden transitions, hence are less 

sxisceptible to scattering effects by several orders of magnitude. The resonance 


f ' 

f 

^ and triplet lines are grouped closely together in wavelength and consequently, 

can be observed with jtdniinal time ca: wavelengtli dependent effects. Because the 
relative intensities of the forbidden transitions are density sensitive, we 
^ have adopted the triplet to singlet intensity ratio (defined by Gabriel and 

I Jordan (1969) as the G ratio) as an indicator of resonance scattering effects , 

I_(ls^ h - ls2p ^P) + I..(ls^ - ls2s ^S) 

G = 21 

Ij^ds S - ls2p "T) 

I 

As the sum of I_ and I_ is constant with electron density for an equilibrium 
X r 

^ plasma, any variation in G will indicate a relative weakening or enhanoenent 

of the resonance line with respect to the intersystem and foi±)idden lines. 

( For exaiiple, if G were observed to be larger than the expected value then this 

would be an indication of a relative weakening of the resonance line. 

i 

Ihe G ratio shows only a weak dependence on electron tenperature. This 
has been studied by Mewe & Schijever (197S a,b) and more recently, by Acton and 
Brown (1978) . Acton & Brown take into account in their calculations the contribution 
from unresolved satellite lines and all excitations to the n=4 levels. Using 
this theoretical value of the ratio (Qp) and oonparing it to an observed value 
(G^) provides a measure of the total change in the resonance line 
intensity, AR, required to produce this effect. If the unscattered resonance 
^ line intensity is then the observed intensity (R^) is given by:- 

R^ = + AR ' (2) 
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It follows that if the forbidden and intersystan lines remain unaffected by 
the scattering prooess then AR is given by:- 



(3) 


A positive value of aR would indicate a net increase in the resonance line 
intensity, where as a negative value would indicate an overall loss of radiation 
due to scattering out of the line of sight. It is not possible to directly 
e>$)ress the optical depth of the region in terms of AR as it is a result of 
several counteracting factors. There will be a loss term due to absorption 
along the line of si^t. However, this can be counteracted by contributions 
to the radiation field from any sources of emission, including the ambient 
corona plasma and nearby active region emission not included in the field 
of 

In the work presented here, we furnish conclusive evidence for significant 
variations in the G ratios of He like ions. In the next section this data is inter- 
preted in terms of resonance scattering following the work of Acton & Brown (1978). In 
the subsequent section a simple coronal model is developed taking into account 
the geometry dictated by the angular redistribution function introduced by 
resonance scattering. This model is applied to the data for ccrtparison. In 
the last section, we discxoss seme possibilities of lasrlng this technique as a 
diagnostic tool in the interpretation of coronal plasmas. 


2. The Data 


In 1973 the Lockheed Palo Alto Research Laboratory flew two sounding rocket 
ejqaeriments that carried an array of six collimated Bragg Crystal Spectrometers. 
Several of these spectrometer systems nonitored the emission frcxn heliiit>-like 
ions in the soft X-ray part of the solar spectrum. The payloads had tlie capa- 
bility to change their pointing positions under direct ej?>eriinenter control. 

Henoe, a series of raster type scans were executed over areas of active regions 
and quiet oorona. 

The first flight, Aerobee 26.002, was on 11th June, 1973, at 
18:30 UT. It took data at a total of 20 pointing positions over the disk. Ihe 
three primary regions observed were: a region of quiet oorona near the disk 

center, McMath 12379 and MCMath 12387. McMath 12379 was a fairly quiescent 
active region near the disk caiter at the time of launch. A total of nine 
observations were rcacte in and around this region. McMath 12387 had just appeared 
over the eastern linb and seenved to be a more oonplex region. A total of eight 
observations were nade of this region. Three spectrometers on this flight were 
dedicated to monitoring He- like ions. The first was an instrument with a wide 
field of view (150" FWHM) , vhich observed the 0 VII system at 21. sS, 21. sS and 
22. iS. A second instrument monitored this system with higher spatial resolution 
(35" FWHM) . The other detector observed the Ne IX system at 13.448, 13.55A, 
and 13.7oSatthe higher spatial resolution (35" IVJHM) . 

The second flight, Aerobee 26.003 took place on 19th Decenber, 1973, at 
22:07 UT. It took data from eleven stable pointing positions in ajid around Mcr'*ath 
12664. This was a very acti-.;e region prior to the fli^t, producing a series of 
flares and subflares. Again three instruments observed He-like ions, but all 
at a spatial resolution of 35" (FWHM). As for 26. 00 2, the 0 VII and Ne IX ions 


were nonitored, the additional ion on this fli^t was XI at 9.178, 9.248 
and 9. 328 . A nore detailed descripticai of the payloads is given by Strong 
and Acton (1979) and a detailed description of the data reduction techniques is 

given by Strong (Thesis, 1978) . Four sarple spectra are shown in Figure 1. 

Using line intensities derived from these data, the G ratios were 
calculated and are listed in Table. 1. The uncertainties, sha^m 
in parenthesis, are cterived from the unoertainties in the individual line 
intensities obtained in the data reduction processes. They probably tend to be 
pessimistic. This is especially true in the case of the G ratios where we are 
calculating the ratios of lines observed by the same detection system and at 
similar wavelengths. Henoe, a ratio of such lines will eliminate absolute . 
calibration errors and reduce to a minimum any wavelength dependent effects. 

The mean values, vhich are mainly vised in the subsequent calculations, are 
derived from the sum of the individual line intensities over a given region. 

The 0 VII results show significant variations from one region to another. 
The wide field of vies? data demonstrate the consistency of the ratios within a 
region due to the integrating effect of a large field of view. The 35" 
field of view 0 VII observations show larger fluctuations, vhich are probably 
caused by a combination of reduoed statistical reliability of the line 
intensities and genuine variations due to changes in the cxjronal structures 
within the various fields of view. The quiet corona has a Icwer value of G 
than any of the active regions. There is an overall indication of an increase 
in the G ratio with a decrease in the size of the field of view. This effec± 
is only jvist within the resolution of the data so caution must be exercised 
in eissigning much significance to it but a possible explanation is discussed 
later. The Ne DC results do not show any significant variations outside the 
experimental uncertainties. 
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Table 2 sumnarizes the average properties for each region for the 0 VII 
results. The effective tenperatures were calculated fran the relative inten- 
sities of the brightest lines observed. The calculations assumed a state of 
(^namic equilibrium existed in the emitting plasma. Although these can only 
be considered to be approximate values (i 0.3 10^ °K) the insensitivity of the 
G ratio to tenperature renders this uncertainty in tenperature to be 
uninpoirtant. The electron tenperature of the quiet corona was found to 
to be 1.9 10^ ®K. The active region tenperatiures are topical of those calculated 
for similar regions from other X-ray data. Using these effective tenperatures 
the ejpected G ratios were calculated from Acton and Brown (1978) . The agree- 
ment is very poor for all the results, none of them falling within the un- 
certainty of observation. 

Alternatively, if we assume that the observed G ratios were caused purely 
by tenperature effects, it is possible to use the ratio to give an indication 
of the tenperatures of the regions. The results of this approad'i are also 
listed in Table 2. This gives the unrealistic result that the quiet corona is 
at a temperature of 6.10^ ®K, and is considerably hotter than active regions. 

Both Meve and Schrijver (1978a, b) and Anton and Brown (1978) discuss non- 
equa librium plasnvis. These can cause changes in the G ratio. However, it would 
seem unlikely that all tliree active regions were in a state of reoarbination 
(especially Mdlath 12379 and 12387) . The quiet corona would have to be in an 
ionizing state to account for the observations. The 0 VIII to 0 VII line ratios 
are consistent with equilibrium values for electron torpratures of 3.5 (i 0.5) 

10^ °K for each acti''.e region observation. If we were to use the sum of the 
triplet lines to derive the tenperature instead of the observed resonaioe 
intensity then this tenperature would be lower. There is, consequently, no 
observational justification for assuming non-equilibrium conditions. 


>-6 


We suggest, in the absence of a viable alternative, that the observed variations are 
due to perturbations of the observed resonance line intensity due to scattering 
of X-rays in the corona. We will subsequently investigate this contention, and 
see if it is viable ^or a set of reasonable coronal conditions. Using (3) it 
is possible to calculate the change in the resonance line intensity required to 
produce the observed effects. In order to increase the G ratio to the levels 
measured in the active regions, the resonance line intensity has to have been 
reduced by between 10-42%. Perhaps even more surprising is that the resonance line 
has to undergo an enhancement of 37% from its theoretical value to reproduce 
the observations for the quiet corona. Ihis is, potentially, a 
difficult result to reproduce from a model, as the region is sipposedly 
of low density, and consequently less subject to scattering effects. 

In order to investigate resonance scattering as a viable explanation of these 
effects, we choose to model the quiet corona case for several reasons. Firstly, 
as nie.ntioned above, it is ^e most difficult region in which to get 
a measurable effect. Secondly, it is a fairly sinple problem geometrically, 
at least in oonparison with the complex loop structures associated with acti\e 
regions. Also, it is a region v^iich is likely to be in equilibrium, stable, and 
could provide the simplest tests of the effect for future observations. 

Results from other spectroTiieters flcvm in recent years 
also indicate variations in the G ratio, these are limited in Table 3. 

These shoiv' variations in the 0 VII G ratio similar to those of 26.002 and 
26.003. The Ne IX G ratios seemtoshav some significant variation when oonpared to the 
quoted errors. A quite interesting result is that for Mg XT which ehcws a 
very large change for McKath 12624. This region was associated with 
a flare at the time of observation, hence the low G value is possibly associated 
with an ionizing plasma rather than a significant scattering effect. 
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3. 


The Model 


We shall assume the oomia to be represented by a spherical shell of inner 
radius and outer radius R 2 with uniform tenperature and density throu^out. 

Our coordinate system, viuch is shown in Figure 2(a) , plaoes the center of the 
sun at (0,0,0) and the observer at the Earth on the +Z axis. An X-ray point 
source, S, is located at arbitrary position (Xg/Vg/Zg) within the shell. The emission 

from, the surface in the direction Z fran a column C with coordinates (>’.,>y ) 

as seen by the observer will be calculated. The coronal plane (X,Y) at 2=0 

as seen by the observer is shown in Figure 2(b) . Since there is rotational 

symmetry in this plane, we arbitrarily place the point source on the +Y axis. 

There is also synretry about the Y axis so that we need only consider positive 

values of x . 

c 

The radiative intensity fran column C consists of three contributions: 

(1) themal erassion from the quiet corona within the column, (2) radiation from the 
source, S, that is absorbed within C and is scattered in the 2 direction and (3) 
radiation emitte from the surrounding quiet corona, absorbed within C, and scattered 
in the 2 direction. As the optical depth of the column approaches zero so will the 
scattered contributions. 

The formulation of the radiative transfer equation for X-rays in the solar 
corona has been oonsicfered by Acton (1978) and we shall make use of his results. 

We shall consider the line absorption profile (f>^ to be (Saussian, i.e.: 



(4) 


where v= (v-Vq)/AVj^ , 

and are the line frequency and wavelength, 

T is the tenperature and M is the ion mass. 



(5) 

is the Doppler half width. 
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is nomalized such that 

V 

00 

/ ^ dv = 1 

0 V 

The line center absorption coefficient for resonance lines is given by:- 


i n^ 


K = 


irc 


Av, 


(7) 


vhere e and m are the electron charge and nass, c is the speed of light, f the 
absorption oscillator strength and n^ is the ion ground state nunber density. 
The line center optical depth, is where a is the path length. It 

is also convenient to define a mean absorption coefficient Ic by:- 


K 

V 




( 8 ) 


The general solution of the transport equation for the- specific intensity 


is:- 


I-. (2) 


2 -Ty(2,z') 

S^( 2 ') e (z‘) dz' 


(9) 


vhere z is the depth coordinate, the absorption coefficient, the source 
function and the optical depth is defined 


T^{2,2') = / (Z' ) dz' . 

2* 


(10) 


The source functiai can be divided into thermal and scattered onnission parts 

s 


c + e 

C - V V 
"v “ 4 “k 


( 11 ) 
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where and e® are the thermal and scattering anission coefficients, respectively. 
Assiraing oonplete frequency redistribution for the scattered photons the solution 
to the transport equation along the 2 axis becomes (Acton, 1978) .* 




TT 2tT 

/ I 


0 


V-(£) sine 




(r.,e)d-:d?d. ' 


+ (n, e) dv' ) 



(z) 


K dz 

V 


( 12 ) 


v^ere V (6) is the angular redistribution function ipon scattering and^, 
and , are the incident intensities at the scattering volume in colurm C from 
the quiet corona and point source. 6 is the polar angle between photons incident 
ipon column C and the 2 axis, v»hile a is the azimuthal angle between the incident 
photon and X-axis. These three terms represent the components discussed abo\’e. 
FOr the type of atomic transitions considered here (j=0, Aj=l) the redistribution 
pt)ase funch:i.on is identical to that for Rayleigh scattering, i.e. , 


1 ^( 0 ) = 3/4 (1 + cos 6). 


(13) 


The function o(y, di^ to emission fran the quiet corona surrounding C 
is given by:- 


/ (’'i,6) = ^ , (n,e,t) (n,9,t) exp[-^/ K^’,(t)dt]dt (14) 


where t is distance neasured along the incident photon direction from 0 at 
the volume element of C under cjonsideration to the edge of the corona, d^. 

With uniform cxjnditions in the corona^ for a volume element at distance z 

2 

* (Acton's Eq. 30 is in error by a faertor of 1/4 tt in the second berm. Also all his 
preceding equations for e® should be multiplied by 4ti) . 


vilong the z axis reduces co: 


(e»2) 


c,,i . 

(1 - e ^ 

4tt< 



( 15 ) 


Due to the synmetiy of the problem, a is not a function of n. Since we anticipate 

the optical thidoiesses to be encountered in the corona at resonance line 

centers generally to be less than 1 it is convenient to ejqjand e”^in a power series 

to be able to integrate Ea. 12 in closed form. The functions f , f , , and f-, 

O i dL 

which will appear ir the subsequent equations, represent the series expansion in 

2 

oonbination with <f^and({^ . Performing the n and v' integrations now gives for 
the quiet corona scattering term in Eq. 12. 



R / I / '^'{e) sin e d 
1 e=o ® 


( 6 ) 


f^ (0) d.6 3 e ^ dz 


(16) 


where e = e, . . 

V V 

For coraparison with the observations we require the frequency integrated 
_ " 2) 

emitted intensity I = / I dv. Ejqpanding e in a pcver series for 

-<^(^2 “ 2 ) 

<^(^2 “ 2 ) <1 and integrating over '•> we find we can replace i}>^ e 
in Eq. 16 integrated over v by ^2^^^ * 

The total amount of energy absorbed by any particular coronal volxme 
element from the rest of the quiet corona is identical for coronal elerents 
at the same radius from the center. Henoe, for purposes of calculation of 
scattering from the quiet corona emission, we place tite column C at the 
X = 0, Y = 0 position. When the incident photon directic® is tangent to the 
sphere R^, the angle 6, which we call 6^, is given by: 
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( 17 ) 


= circsin (R,/z) . 
m 1 

Ftor e < 6 the distance d is the distaiioe to the sphere vhile for 6 > 6_ , 

in C 46 

d is the distance to the sphere R, . The anount of quiet corona radiation 
c J. 

scattered by a volume element at radius r into an angle 6 from the outward 

9 

radius vector is: 



(1C) 


vhere 

dj^ = 2 cos e +y sin^ e 

= -z cos e + “ 2^ sin^ 6 (20) 

In order to obtain the quiet corona scattered intensity we integrate the 
ccaitribution function from (16) along C with the appropriate r, 6^ and 6^^ for 
each volume element. 

For a point source of intensity I ' located at x , y , z the function 

s s s s 

can be written 


■ ^s' *v' ■ • 'v' >^3 ' 


( 21 ) 


vhere d^ is the distance between a point on column C and the point source. 
Letting !_' = I_ , and integrating over v'we find 

5 9 V 


= ig - 

0 K 


After integrating over outgoing frequency, 
the point souroe scattering tern in the transport equation new beccnes, 


( 2 „ - Z)‘ 


-sB.) ■ -S — J '=> *2 

2d d_ 


vdiere = I /47rd^, and is the total souroe strength, 

s so g so 


Ihe upper limit of integration is the intersection of the column with the 
outer shell and the lower limit 2 g is the intersection with the inner shell. 
2 ^ and 2 g are given by: 


o D 2 2 

= R~ - r 


2 „ 2 2 

"b ' -■'c 

2 2 2 

where r„ = 

c c c 


For many oonfigurations the entire length of the column C is not illuminated 
the point souroe as the souroe is cutoff by the shadowing of the inner shell. 

In this case, as illustrated in Figure 2(a) the lower limit of integration becones 
defined by the intersection of the line tangent to the shell from the source 
with the column C. Appropriate ejpressions for have been derived and 
applied in the calculation. 
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. r 


A ooirplete solution can be obtained by moving the point source along the 
positive y axis and calculating the emitted intensity for all x > 0. 


Itie expression for the integrated emitted intensity now becomes: 


— ^A ^ i ^ 

4ir I e f ^2 2“ •’ ie'*" *^'s ^^2 ^o 

% % 


(26) 


IT , 31 2a [d^ + (z - 2)^] 

+ / »3 Sine de I dz * J -j f, f, <te 

m '*‘0^ 


"2 


S' 4tt 


I (I - ii) -Jig -Jig I 


viiere the first term is the thermal emission emitted directly frcxn C(x ,y ) , the 

o c 

second tem is quiet corona radiation scattered by C ejnd the third term is radia- 
tion from the point source scattered by C. See Figure 2(a) . 
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IV. Numerical Results for 0 VII 


As discussed above the observed X-ray signal f ion our model- corona consists of three 

components! (1) direct thermal emission from the point source and the colurm 

of oorona under observation; (2) scattered radiation from the point source 

and (3) scattered radiation from the quiet corona. The major physical 

dependencies are as follovv's: Conponent (1) is a function of e; corrponent 

(2) is a function of and ocarponent (3) is a function of e and The X-ray 

2 

signal from a one or. column of corona is shown in Figure 3 as a function of 
radios frcm the disc center. For the quiet oorona we assume T^ - 1.9 10^ °K and 

p 

n^ = 2.4 10 . The increase in coirponents (1) and (3) frcm center to limb basically 

results fron the increase in length of the column C. Qmponent (3) has a 
relatively larger increase between center and liitb as the average scattering 
angle becomes more favorable along with the increased column length. At 
the limb, r = R^, the values essentially double as we "see" through a double 
length of oorona. Ccirponent (2) is illustrated by point sources located on 

2k 

the R, surface at various distances r„ from disc center (z-axis) with I ^ = 10 7 
is so 

2 

Ihe intensity falls off as 1/r from tlie source. 

Figure 4 illustrates the scattered radiation from el ments at different 
heights, R, in the quiet oorona as a function of the angle, 6, between the 
radius vector to the elanent and the line of sight. These variations are 
equivalent to a change in the observed scattered signal (in arbitrary units) as 
the elanent moves with respect to the central meridian (6=0) . When the scattering 
element is located at the base of the corona (R=Rj^) only a relatively small pro- 
portion of the atmospheric shell is above the local horizon and able to contribute 
to the radiation field at the elanent, 4 (a) . As R increases, the volume of oorona 
visible from the scattering element increases and therefore so does the signal in 
general. The maximum ''visible'* volume in close proidmitv to the scattering 
element occures for R v (Rj^ + R2)/2 vdiich is vhere the maximum signal occurs, 4(b). 

Ihe rraximum angular scatter for each R occurs in the direction in which the largest 
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fraction of photons are incident. For R = this angle is slightly less than 
ti/ 2, vJiile for R » R 2 the angle is slightly greater than n/2. As we v^ould 
ea^ct in light of the j^ase function, the ooiresponding intensity miniina are 
found n/2 away from the rraxixna. 

Fig. 5 illustrates the disc center signal variations with electron density 

■7 8 

for electron densities between 8 10' and 8 10 . The temperature is held 

fixed at 1.9 10® °K with ^2 ~ ^1 * scattering column is located 

at a point one arc minute from the disk center source region. The point source 
scattering is dependent only upon the line absorption coefficient which varies 
linearly with n^. The thenral emission varies as n^ in direct proportion to the 
volume emissicxi coefficient. The quiet corona scattered signal is proportional to 
both the thermal anission ooefficiait and the line ^sorption coefficient and 

3 

hence varies eis n . In all cases the curves are not linear due to optical 
e 

depth effects. Ihe observed source G ratio is also shown in Fig. 5. 

-k 

The line absorption coefficient varies as T . Since the scatter of a 
particular point source depends only on the absorption coefficient, this scattered 
signal decreases as the tenperature increases. This effect is shown in Fig. 6 
where we also show the variation of the other oonponents and the G ratio. Ihe 
volume emission coefficient increases with T and hence so does the thermal 
anission. The quiet corona scattered ccrtponent is tiie product of the above 
two dependences. 

We ooirbine the above effects to produce the ocxitour plots shown in Figure 7. 
In order to illustrate the relative inportanoe of the scattered ocnponents in the 
observed 0 VII resonance line signal we define a quantity, A, such that. 


A = 



(27) 
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v«her€ and are the total predicted and thenral fluxes in the resonanoe line 
respectively. This particular form of A is chosen as it readily demonstrates the 
variation of the scattered signal as a function of airbient thermal signal, and 
its sign indicates vhether scattering into or out of the line of sight is domi- 
nent at any given location. We show four distinct cases for different source 
positions, from the disk center (7a) to the liirb (7d) . The value of delta 
obtained by observing the source directly is denoted by As, and is primarily 
deteimined by the oolunn length of overlying quiet coronal iraterial. As is 
negative which indicates a net loss of resonanoe line radiation. Ihe contours 
are plotted at suitably spaced values of A in the X-Y plane, cm a scale of a 
unitary solar radius. Ihe limit of the irradiated volume of corona by the source 
is indicated by a broken contour. The choice of appropriate moctel parameters 
is discussed at length in the next section. 

Figure 7a shows the source at the disk center (0,0) . This case has oonplete 
synmetr^' in X and Y directions. For a soxoroe at the base of the corona as is 
-0.14. Otherwise on the solar disk a is positive indicating a net addition 
of scattered radiation along any givei line of si^t. Iimediately above the 
limb, vhere an extended column of plasma is present along the line of sight, 

A becomes negative, but as we look higher in tie atmosphere A tends to zero as 
the CO limn becomes shorter. As the distance of the source from the disk center |] 
increases so does the depth of the overlying anbient corona, hence As beoomes 
mere negative until it reaches a limiting value of -0.43 at the linb. The 
syinietry is maintained almg the X directiem but not along the Y direction as 
the columns available to scatter equidistant frem the souroe in the ordinate 
are different. They are longer being further frem the disk center due to the curv 


ture of the sun, hence will give rise to a greater scattered signal. As the so; 


approaches the liirb the negative contours above the lint) becxane distorted 
(7c) / until finally when the source is at the liirb (7d) tha effect of the 
opaque inner shell beooines apparent producing a sharp horiaon effect. Eelow 
the source horizon the optical depth effects keep h negative as in the first 
case 7a) where as within the sphere of influence of the source L is primarily 
positive due to radiation being scattered into the line of sight from the source. 
In reality one would not ej^)ect such a distinct boundary for the source horizon 
especially if the source region, assuming it to be an active region, is extended 
and not at the base of the corona. However, these models may be used to de- 
termine the gross characteristics of scattering in the oorcna. 
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Section V» Oonparison of the Scattering Mjdel to Cfcservations 

Itie problem that arises in applying this nodel to the observations cited 

in section II lies in the physical parameters to be adopted to represent the active 

region and quiet oorona without predisposing the result. The quiet coronal 

teiperatures can be derived frar the relative intensity of the X-ray emission lines 

observed by the 26.002 spectrometers. It is 1.9 (±.2) 10^ °K (Strong, 1978). 

From the predicted intensity of the 0 VII resonance line, calculated froti the 

triplet lines and the e:q>ected G for that tenperature, the emission measure 

required to give the observed flux can be estimated. Knowing the size of the 

field of view, it is possible to derive a quantity that is the emission measure 

2 

cdong the line of sight for a column of 1 cm cross-section. This has a value 
26 “5 

of 5.8 10 an . A family of heights and densities can be adopted that would 
reproduce this value. It should be noted that this quantity is particularly 
sensitive to density; halving the density would require a four- fold increase 
in the extent of the atmospheric shell. Hence, the choice of a representative 
height for the shell would restrict the density to a relatively narrow range of 
values. On the basis of the extent of liirb brightening on soft X-ray photo- 
graphs, and the variation of coronal density and tenperature wd.th height (Allen, 

1973), we adopted a value of 1.10^*^ on. This determines the electron density 
8 —3 

to be 2.4 10 on . This also determined both the emission and line center 

absorption coefficient? to be 72 photon/cmVsec and 1 59.10 an ^ respectively. 

49 

This model corresponds to an emission measure of 4.10 for the general corona 

49 

and oonpares favorably with other determinations vhich lie between 3.4- 4.4 10 
(Allen 1973; Acton et al. , 1972) . 
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Aeiobee 26.002 took observations of both active regions and the quiet oorona. 

The quiet coronal G ratio for 0 VII (0.63) is low oonpared to the theoretical 

value derived from the observed tenperature (G = 0.86-.88). Hje observations 

also give the position, taiperature and emission measure of the major region 

near this field of view, fron these data it is possible to reoonstruct an 

32 2 

equivalent point source intensity for NcMath 12379 of 6.0 10 photons/cm /sec/ 

St. The coordinates of the active region, along with those of the pointing posi- 
tion for the quiet oorona observations, were luicwn frcm the aspect solutions of 
the rocket, and confirmed from Ha photographs. The source was put at a height 

9 

of 5 10 an in the atmospheric shell based on the size of the loops structures 
derived from ATM photographs taken in the X-ray and EUV spectral regions of McMath 
12379. 

We have defined a quantity, L, vhich reflects the iinportance of scatter 
at any position in X and Y. The case under consideration above is closely 
represented by Figure 7b. .To ejsplain the observations of the quiet corona made 
by Aerobee 26.002, A would have to have a value of 0.27 (±0.11) . Hds model 
predicts a value at the required position of A = 0.25 using the above model 
parameters. This demonstrates that scattering can be a significant phenomenon 
in the aitbient oorona euid that its magnitude is sufficient to explain observed 
changes in line ratios such as the G ratio. Ihe ability of this simplified 
model to predict the observed change in the resonance line flux and consequently 
the anomalous G ratio for the quiet oorona cbservation is very encouraging, 
especially without any adjustment of the model parameters from tl:ieir independent 
values. Although the scattered signal is density sensitive, any change in 
density would have to be accompanied by a corresponding change in the extent of 
the atmospheric shell to maintain the emission measure constant. Ihese 
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oKXJsing effects will tend to oonpensate for one another making the model less 
sensitive to changes in density. Figure 8 shows the variation in the predicted 
A with a change in the extent of the corona while maintaining the total emission 
measure (hence the observed flux) at a constant value. It can be seen that the 
model is relatively insensitive to such changes in ocmparison to the uncertainty 
on the observed value of A. 

Ihe question remains as to whether scattering effects from the quiet 
corona are responsible for the deviations in the observed active region G ratios . 
Figure 9 shows the results of the calculation of the source region G ratio as a 
function of the radial position of the source. The source is assumed to be at 
the base of the atmosphere. The G ratio rises rapidly as the source nears the 
limb and then falls as the source is moved higher in the atmosphere above the limb. 
Tnese effects are due to the cliange in the length of the column of quiet coronal 
material along the line of sight to the source. We ass'ume tiie source G ratio to be 0.75 
without scattering. Plotted on this graph are the G ratios derived from the large 
field of view data (150" FVJHM) frcm 26.002 for both active regions. VJhile neitiier 
groips of points lie exactly on the curve the general trend is in good agreement. 
Adjusting the curve to fit either data set more closely by changing the quiet 
corona density, for exanple, will result in a poorer fit to the other data. 

The small field of view ratios do not shew such a good correlation but the 
uncertainties and fluctuations are too large to draw any firm conclusions in 
this respect. However, the mean values over each region shav a similar trend to 
that of the large fields. The fact that these observations do not fit is 
probably caused by the assunption of a point source and that it is assumed to 
be optically transparent to the 0 VII radiation. This will not be the case, and 
the emergent radiation from the dense active region loops might give a different 
G ratio for 0 VII that the valxae of 0.75 assumed above due to scattering 
within the locp structure. 
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As noted in Section II, the larger fields of view tended to slww a lcv;er 
value for G. This c»uld be due to either the inclxision of an increasing amount 
of quiet corona emission (with a different G ratio characteristic of its tempera- 
ture) or incli>ding a proportionally larger amount of scattered signal. We 
investigate these effects in Figure 10, vhere they are contrasted. The source is 
placed at the disk center and the field of view observes the source directly. 

As the full width at half maximm (PWHM) of the instrument is increased, the 
resulting G ratio is calculated with and without scattering effects included. 

The model parameters are as above for the quiet corona. The source intensity 
is a variable, and the source is at the base of the cx>rona. 

For the no scattering case the G ratio is seen to increase as a function 
of the field IV5HM at a rate that is related to the source intensity (solid 
curves) . This is due to the increase in the prc^rtion of the signal from the 
lew temperature quiet corona which has a higher G ratio. With scattering effects 
also included, the G ratio drops as the PWHM increases as is shown by the broken 
curve, again the rate is in proportion to the source bri^tness. Ihe starting 
value of G is determined by the column of overlying quiet corona, for this case 
(R = 0) it would give a value of 0.88 for G. If the source were at the linb 

5 

34 

starting value of G would be 1.35. The source brightness is assumed to be 10 
2 

phts/an /sec. The reuslts are oenpared to the data from Mc3^5ath regions 12379 
and 12387 vhere we have two spatial resolutions available. The overall trend 
of the observations confirm the scattering model over the von scattering model. 
However, this effect is narginal in these data and caution should be ei:ercised 
in interpreting these results. 
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VI. concluding Discussion 


We have presented the observed G ratios of several helixim like ions 
originating from a large nurtber of discrete pointing locations over the solar 
corona. We have noted that the ratios have a large spread of values, and that 
the anomalous variations are too large to be explained in terms of the 
statistical fluctuations in the data or from physical differences in the emitting 
plasma, such as in tenperature or equilibrium balance. We conclude that scattering 
of X-radiation in the corona is responsible for the observed variations. We, 
then, investigate this supposition by developing a siitple coronal model and 
applying it to the observ^ed situations for the 0 VII results. The major 
variations are observed in 0 VII, vhidi woiold be expected as this ion is the 
most susceptible to radiative transfer effects at coronal temperatures. 

The radiative transfer code developed during this work takes account of the 
Raylei^ angular redistribution function and the geometrical factors associated 
with the model. Ihe atmosphere is assumed to be an isothermal homogeneous 
shell. The thickness, density and tenperature of this shell may be changed 
independently. A point source of arbitrary brightness can be placed at any 
coordinate and hei^t within the shell. The relative contributions from the 
thermal enission from the source and quiet corona are calculated as seen ty 
the observer over a matrix of points. These represent the signal that would 
be expected from each elertent if there were no scattering. The scattered signal 
from the source and surrounding ambient corona is also oorputed. This can be 
a positive or negative quantity depending vrfiether there is a net loss or gain 
of resonance line radiation. If there is a net gain of resonance flux, for 
exaitple, then the G ratio would be depressed ocnpared to its theoretical value. 
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We use the itodel to determine the sensitivity of the scattering mechanisn-,, 


and consequently the G ratio, to physical parameters of the plasira. We 
;^ly the model to the problem of a volume of ambient corona illuminated by 
a nea3±>y active region. The tarperature, density and extent of the observed 
region, as well as the source brightness, were derived in(^pendently of the model. 
In this case, the predicted scattered signal fron the model agreed well with 
the reguireonent from the observations. Ihe result was insensitive to the 
detailed form of the model. However, an attenpt to interpret the active region 
observations using this model were not as convincing as the previous case. 

However, the basic trends in the observed G ratios are reproduced by the model. 

We believe tlie discrepenc^' to be due to the inadequacy of the assunption of the source 
region being a point source and that it emits radiation characteristic souly 
of its tenperature, i.e. a G ratio of 0.75. This assunption is equivalent to 
presuming that there is no scattering within the active region itself, which 
will probably not be correct. A stud;' is underway to investigate the problem 
of scattering within active region loops. In this case, the geometry will be 
the dominant factor, but it will help us to understand the observations from 
the hi^er spatially resolved data, vhere individual structures become inportant . 

The phenomenon of resonance scattering with its sensitivity to density, 
iray be a useful diagnostic tool in tlie future. Conclusive evidence for its 
presence in the corona could be furnished by the detection of a polarized 
oonponent in the observed X-ray line intensities. Where -ever there is an 
anisotropy in the radiation field, as would be provided by' a nearby source 
region, then the emergent radiation would be partially polarized. As the crystals 
used in Bragg crystal spectrometers produce a polarizing effect a change in the 
aspect of the instrument with respect to the sun will modulate the observed 
signal. This would be observable given a sufficiently sensitive instrument. 
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In conclusion, we believe that resonance scattering is an inportant 
mechanism in the solar corona, at least for sane X-ray anission lines. Ihis 
should sound a note of caution when using X-ray line intensities to derive 
nodels of the coronal pleisma. These results show that substantial errors can 
be introduced in some sitvitations if the optical depth of the medium and its 
geometry are not taken into account. 
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TABI£ 1. Observed G Ratios of Ke-like Lines in Coronal Spectra . 


Pointing Position 

0 VII 

0 VII 

Ne IX 

Mg XI 

Label 

Aerobee 26.002 

150" FWIM 

(35" FVHM) 

(35" PWHM) 

(35*' FWHMJ 

CTR PRE 

0.62(.05) 

- 

- 

- 

CTR RL 

0.72(.07) 

- 

- 

- 

CTR PST RL 

0.56(.04) 

- 

- 

- 

Mean Quiet Corona 

0.63(.07) 

0.89(.40) 

- 

- 

Move to Mctoth 12379 

0.79(.03) 

- 

- 

- 

AR131-1 

0.85(.01) 

0.65(.15) 

0.95(.17) 

. 

API 31- 2 

0.S3(.01) 

1.13(.20) 

0.75 (.12) 

- 

AR131-3 

0.86(.02) 

1.03(.14) 

C.89(.ll) 

- 

AP131-4 

0.86 (.03) 

0.84(.12) 

1.55(.58) 

- 

AR131-5 

0.77(.05) 

1.54(.51) 

- 

- 

AP131-6 

0.89(.03) 

0.77(.49) 

- 

- 

AR131-7 

0.82(.01) 

0.83(.40) 

1.32(.63) 

- 

AR131-8 

0.85(.03) 

1.27(.31) 

> 0.95(.40) 

• 

Mean Mcr»3ath 12379 

0.84(.04) 

0.96(.15) 

0.98(.16) 

- 

Move to McMath 12387 

0.8K.05) 

- 

- 

- 

AR137-9 

1.06(.02) 

1.39(.24) 

0.97(.18) 


AR137-10 

1.09(.01) 

0.96(.46) 

- 

- 

AR137-11 

1.12(.02) 

1.30 (.21) 

0.85(.24) 

- 

AR137-12 

1.06(.02) 

0.99(.14) 

0.82 (.20) 

- 

AP137-13 

i.07(.01) 

1.68(.75) 

- 


AR137-14 

1.04(.02) 

1.29(.18) 

1.06 (.13) 


l^an McMath 12387 

1.06(.03) 

1.22(.20) 

0.96 (.20) 

- 

Aerd^ee 26.003 

1 

- 

0.99(.08) 

0.80(.07) 

1.34 (.20) 

2 

- 

0.89 (.07) 

0.82(.08) 

1.02(.16) 

3 

- 

0.68(.06) 

0.98(.ll) 

1.02 (.19) 

4 

- 

0.66(.30) 

1.24(.20) 

- 

5 

- 

1.16(.21) 

0.97(.15) 

- 

6 

- 

0.89(.08) 

0.89(.ll) 

1.34(.34) 

7 

- 

1.14(.21) 

0.70(.06) 

1.19(.20) 

8 

- 

1.19(.15) 

0.75(.ll) 

0.99(.33) 

9 

- 

0.73(.06) 

1.11 (.28) 

- 

10 

- 

1.06 (.24) 

0.86(.14) 

- 

11 

- 

0.64(.30) 

- 

0.96 (.29) 

Ptean Mctoth 12664 

- 

0.94(.12) 

0.84 (.109) 

1.14 (.20) 
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Table 3 


Results from Other Measurements of the G Ratio for Various He-Like Ions 



rov 

( "RnWI) 

Ion 

Active Region 

0 VII 

Ne I.X 

Mg XI 

Source 

McMath 11060 

240 

1.27 

0.89 

0.93 

(b) 

McMath 11619 

180 

0.87(.13) 

0.75(.04) 

0.77(.16) 

(b) 

McMath 11621 

180 

1.13(.05) 

0.88(.03) 

1.04(.05) 

(a) 

McMath 12624 

180 

1.25(.14) 

0.95(.01) 

0.69{.03) 

(a) 

McMath 12628 

180 

0.92(.09) 

0.9K.03) 

0.83(.04) 

(a) 

Quiet Corona 

169* 

1.16 

1 

1 

- 

(c) 

1 

McMath 12336 

45 

1.2K.3) 

0.93 
1 

- 

(d) 

; 


* 1-D oollimation. 

(a) J.P^’e, Private Conmuni cation. 

(b) J. Parkinson; 1972, 1975, Thesis. 

(c) McKenzie, et al.: 1977. 

(d) Davis et al., 1978. 
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FIGURE CAPriaiS 

Figure 1: Sanple 0 VII Spectra fron each of the observed regions. The sum 

of tile forbidden and intersystan lines have been us^ as a scaling 
factor in plotting the histograms, henoe the resonance lines can 
be octTpared directly. 

Figure 2a: The coronal model (in the y-z plane) adopted in the calculations. 

The corona is represented by a isothermal homogeneous tenuous shell 
of plasma overlaying an opaque sphere of radius, I^. The source can 
be placed at an ai±;itrary position within the shell. 

2b; The grid (in the x-y plane) represents the matrix of points at 

vhich the emergent flux is oorputed in any given X-ray emission line. 

Figure 3: A ocnparison of the various oonponents that contribute to emergent 

flux obserwd from any position on the sun. The source scatter is 
calculated "is a function of its radial position on the disk. This 
quantity may be scaled according to the source intensity. The thermal 
emission from the anbient corona and its corresponding scattered 
ocHiponent are also shcv;n without scattering effects. 

Figure 4: The ocrrbined result of the angular redistribution of the photons 

with the horizon effects, as a function of altitude in the 

atmosphere for an elemental volume. The abissa is the angle between the 
radius vector to the element sjid the observers line of si^t (z axis) . 

The ordinate is the relative scattered intensity wei^ted by the 
relative brightness of the radiation field at the element. 
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Figure 5: 


Figure 6: 


Figure 7: 


Figure 8: 


Figure 9: 


Figure 10: 


Ihe variation of tlie source scatter, coronal therml emission, and 
coronal scatter as a function of tl'e electron density at a point 1 

arc minute away from the souroe. The source is at the apparent disk 

32 

center, at the base of the corona>and has an intensity of 6.10 phts/ 
2 

cm /sec/st. Ihe resulting change in the G ratio is also shovTi. 

The variation of the source scatter, coronal thental enission and 

coronal scatter as a function of the airbient electron tenperature of 

the corona , at a position 1 arc minute away from a source at the 

apparent disk center. The resulting G ratio is also calculated. 

Contour plots of A for four different souroe positions, (a) Souroe 

at the apparent disk center, (b) Source region 400 arc seconds frcm 

the disk center, (c) Souroe region 800 arc seconds from disk center 

(d) Source at the limb. Ihe souroe is at the base of the corona 

32 2 

and at an intensity of 6.0.10 phots/cm /sec/st. Note the variation 
of the souroe A as the line of sight column density' increases. 

By maintaining a constant emission measure, a series of density and 
height model fit the observed fluxes. Tlve resulting variation of A is 
calculated for a series of models around the adopted value (+) . 

A comparison of the observed and calculated G ratios for the large 
field of view (150" R^HM) data as a function of its radial dis- 
placement. Ihe vinscattered G ratio is assumed to be 0.75. 

A comparison of tlie observed variation in the G ratio v^dth the size 

of tlie field of view (FV2iM) with (a) tlie G ratio for different 

2 

souroe brightnesses (phts/cm /sec) without any scattering, and (b) 

34 2 

the G ratio variation for a source brightness of 10 phts/cm /sec 
with scattering effects included. 
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Appendix 4 


Contents 

Summary of computer programs developed during course of this investigation. 

by: Keith T. Strong 

John L. Kulander 

Key: to Program ID code 

P : Prediction 

D : Data Manipulation 

A : Analysis 

T : Theoretical Studies 

M : Mathematics 
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1. EMITS 


P 


KTS 


Processes any given transition to calculate its temperature and emission function 
It also plots the resulting emission function. It can process any number of 
lines from a line list file between given wavelength limits at one time. 


2. CJBAL - D KTS 

Ionization Balance data file. Elements:- C, N, 0, Ne, Mg, Si, S, Fe 

(IX - XX\’II), Ni (X - XXIX), Ca (IX - XXI). 

3. GLEEM - P KTS 

Faster version of EMITS (1). Options to plot multiple curves of either 

the emission or temperature functions. Utilizes Spline routines to smooth 
plots rather than straight line interpretation of EMITS. 

4. CLONE - A + P KTS 

Statistical analysis of 0 VIII : Ne IX line ratios to determine the 

relative elemental abundance. 

5. POINn and POINT - T JK (KTS) 

Program to study the radiative transfer of X-rays in the solar corona 

for various geometries. It calculates the thermal and scattered signal from 
an elemental column situated in an emitting spherical shell of the n>»;"nsions 
of the corona. 

6. SPLYN - M - KTS 

Interpolation and integration routine utilizing spline techniques. 
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7. DEMON - P 


KTS 


Predicts line intensities given the transition and wavelength by 
-•♦searching a data file. However, will accomodate 

; differential emission measure models in the range of Log(T) ■ 5.5 - 7.0. 

8. MODEL - A KTS 

Takes any set of line intensities and attempts to produce a differential 

emission measure model using a method described by Withbroe (1975) . It also 
plots the resulting model. The program is interactive so the analysis may 
be interogated at any point to investigate any problems. 

9. LOPl - P JK 

Predicts emission observed from a 3-D loop structure of constant cros- 

section as seen from any orientation. Development routine for studying 
resonance scattering of X-ray in a coronal loop. 

10. NRLTP - A KTS 

Takes NRL digitized images and convolves them with the collimator res- 
ponse of the spectrometers deriving an equivalent intensity in the chosen 

EUV lines. It takes the film densities on the tape and applies the instrument 
calibration at any given position in the image. 

11. IMAGE - M KTS 

Takes any function plotted in 2-D and represent it as 3-D image using 

Euler angles. Develope routine for X-ray resonance scattering. 
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